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ABSTRACT 

In a search for evidence of the short wavelength increment in the Sunyaev-Zel'dovich (SZ) 
effect, we have analyzed archival galaxy cluster data from the Sub-millimetre Common User 
Bolometer Array (SCUBA) on the James Clerk Maxwell Telescope, resulting in the most com- 
plete pointed survey of clusters at 850 /im to date. SCUBA's 850 /im passband overlaps the 
peak of the SZ increment. The sample consists of 44 galaxy clusters in the range < z < 1.3. 
Maps of each of the clusters have been made and sources have been extracted; as an ancil- 
lary product we generate the most thorough galaxy cluster point source list yet from SCUBA. 
Seventeen of these clusters are free of obvious AGN and have data deep enough to provide 
interesting measurements of the expected SZ signal. Specialized analysis techniques are em- 
ployed to extract the SZ effect signal from these SCUBA data, including using SCUBA's 
short wavelength band as an atmospheric monitor and fitting the long wavelength channel to 
a model of the spatial distribution of each cluster's SZ effect. By explicitly excising the exact 
cluster centre from our analysis we demonstrate that emission from galaxies within the clus- 
ter does not contaminate our measurement. The SZ amplitudes from our measurements are 
consistently higher than the amplitudes inferred from low frequency measurements of the SZ 
decrement. 

Key words: cosmology: cosmic microwave background - cosmology: observations - galax- 
ies: clusters: general - submillimetre 



1 INTRODUCTION 

The Sunyaev-Zel'dovich (SZ) effect jSunvaev & Zerdovichll97^ 
is a distortion of the cosmic microwave background (CMB) spec- 
trum seen through a reservoir of hot plasma such as is found in 
clusters of galaxies. Measurements of the SZ effect in clusters 
can be used to determine the physical conditions within the intra- 
cluster medium. The SZ amplitude, in combinati on with other data , 
can act as an important probe o f cosmology l lBirkinshawl 1 1 9991 
ICarlstrom. Holdei- & Ree"s3l2002h . 

There are several motivations for measuring the SZ effect. 
Blank field searches for clusters at all redshifts allow one to con- 
strain the cluster redshift distribution N{z). This is anticipated to 
be especially important in understanding the epoch of early clus- 
ter formation since the SZ effect remains bright at large redshifts. 
These studies provide information about the massive end of the dis- 
tribution of collapsed objects in the Universe. Surveys of catalogues 
of known clusters allow one to measure large scale flows via the ki- 
netic SZ effect. This effect is small in any single cluster so larger 
surveys are called for. Unlike the intrinsic SZ effect, the kinetic 
effect is spectrally broad and can be either positive or negative, so 
measurement over a wide frequency range is desirable. In principle, 
a combination of the data presented here and SZ measurements at 



other frequencies provide a measurement of the kinetic SZ effect, 
although in practice this is quite difficult, even with measurements 
at a number of wavelengths (Benson et al, 2004). Detailed targeted 
observations of a single cluster allow the study of physical details 
of the baryon distribution and dynamics in the cluster, especially 
in conjunction with good X-ray maps. These motiva t ions a re all 
discussed in, for example. ICarlstrom. Holder & Reesj | |2002|) . 

SZ decrement measurements at frequencies of tens of GHz 
now seem almost routine. However, detection of the increment has 
proven to be considerably more difficult. The data discussed here 
can be used to address the physical parameters of clusters on an 
individual basis, and are especial ly powerful when com bined with 
data at several other wavelengths l lLaRoque et alj2006l) . Such com- 
binations exploit the characteristic spectral shape of the SZ effect 
to increase the confidence of the overall measurement, and perhaps 
to constrain smaller spectral modifications, i.e. from the kinetic 
SZ effect. Of practical interest, these archival SZ effect measure- 
ments provide insight into the practical difficulties associated with 
measurement of low surface brightness, extended emission in the 
mm/sub-mm regime. 

Ground based sub-millimetre (sub-mm) observations of the 
SZ increment are possible because the peak of the SZ spectral dis- 
tortion occurs near 370 GHz (810 /im) in an atmospheric window. 
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However, while decrements in emission are rare astrophysically 
and can be ascribed to the SZ effect with little ambiguity, the sub- 
mm emission from clusters might not be due to the SZ effect alone. 
Three possible contributions to the flux are: 

(i) the SZ effect itself, 

(ii) gravitational lensing of high redshift background sub-mm 
sources, and 

(iii) dust-obscured starburst galaxies within the cluster itself. 

Despite being foregrounds from the point of view of SZ effect mea- 
surement, items (ii) and (iii) are interesting topics in their own 
right. Although gravitationally lensed background sources offer 
little information about the lensing cluster itself (besides a con- 
straint on the mass surface density in the lens), they allow impor- 
tant measurements of high redshift sub-mm sources which would 
otherwise be too dim to resolve. Virtually all massive galaxy clus- 
ters which have been studied host at least one strongly lensed 
sub-mm source. Several studies have made use of galaxy clus- 
ters as natur al telescopes to pr o be the high redshift sub-mm pop- 
ulation (e.g. ISmail et alj |l99 7l IChapman et al.ll2002l. lBest|[200a. 
ICowie. Barger & Kneibll2002llSmari et alj|2002h . Sub-mm sources 
within a galaxy cluster are rare and, when present, ten d to be as- 
socia t ed with an AGN in th e central cluster galaxy (e.g. lEdge et al.l 
ll999llChapman et al.ll2002h . There are weaker but more problem- 
atic effects due to increased inc idence of dusty star-forming g alax- 
ies at higher redshifts (see e.g. lBeslll2002l . lWebb et al.ll2005f) . per- 
haps related to the Butcher-Oemler type phenom enon seen outsid e 
the cores of some clusters in mid-IR studies (e.g. lDuc et alj|2002h . 
As discussed later, both lensed sources and emission from galaxies 
in the clusters themselves are problematic foregrounds and provide 
some of the fundamental limitations to measuring the SZ incre- 
ment. 

SCUBA on the JCMT, which began common-user operations 
in 1997, has the most extensive archive of sub-mm measurements 
of galaxy clusters available. We have attempted to measure the SZ 
effect increment in these clusters using this data set following th e 
techniques we developed in an earlier paper jZemcov et al.ll2003h . 
Previous analyses of these data have been performed by several 
different groups, and so have been quite heterogeneous in terms 
of the analysis methods employed. In contrast, this paper presents 
the results of using a consistent analysis pipeline for these galaxy 
cluster data. In particular we take great care to remove the effects 
of common mode atmospheric emission in a way which does not 
also cancel the SZ emission. By addressing both source identifica- 
tion and the SZ effect, and by handling all of the data consistently, 
we have produced a uniformly selected list of candidate sources in 
these cluster fields. 

In this paper, we first present the data selection procedure 
and data analysis leading to the individual cluster maps. This is 
followed by a presentation of the SZ effect analysis method and 
results. The point source list and SZ effect amplitude results are 
then discussed. Finally, we discuss some of the difficulties asso- 
ciated with SZ increment measurement in the presence of strong 
atmospheric emission and lensed background sources, and lessons 
learned for future experiments. This paper also contains 3 appen- 
dices: the first presents listings of integration times and weather 
conditions for each cluster field; the second presents a complete 
list of point sources; and the third is a discussion of the characteris- 
tics of each cluster, with a list of references for relevant prior work 
on each cluster field. 



2 DATA COLLECTION & ANALYSIS 

This work has made use of both JCMT archival data and new 
JCMT observations; both are discussed in detail below, followed 
by a description of our analysis methods, which differ from stan- 
d ard SCUBA analysis t ools. These methods follow those described 
in lZemcov et al. to which the reader is directed for more de- 

tail. Essentially, this analysis differs from the canonical SCUBA 
pipeline in two major ways. Firstly, SCUBA's short wavelength 
channel is used to subtract atmospheric contamination from the 
long wavelength (science) channel. Secondly, we fit the set of dou- 
ble difference measurements to a model of the SZ effect rather than 
attempting to image it directly. The motivation for these steps is to 
account for loss of flux due to SCUBA's spatial differencing and 
standard 'sky-removal' techniques. Before attempting to measure 
the SZ signal we identify and remove contaminating point sources. 



2.1 Archival data selection 

A candidate cluster list has been compiled as follows. Firstly, the 
list of proposal titles and abstracts submitted to the JCMT over 
the period 1997 to 2005 was searched for projects likely to target 
galaxy clusters and the target lists of these proposals was examined. 
This forms the core of our target list. In addition, the literature was 
searched for published SCUBA observations of galaxy clusters; the 
list of clusters found in this way completely overlaps with the pro- 
posal list. The final check is to find the complete target list for any 
project that targeted even a single galaxy cluster, which should find 
essentially all unpublished observations. Only unpublished obser- 
vations of cluster fields from projects whose proposed observation 
programme was unrelated to galaxy cluster measurements will be 
missed using these search criteria. 

The on-line JCMT archive^ were searched for data collected 
within a 10 arcmin box centered on the candidate SZ clusters. 
Almost all SCUBA g alaxy cluster data taken in unpolarized 'jig- 
gle mapping' mode ( [Holland et al] 1 19991 [Archibald elal] l2002l 
.Zemcovetal. 2()o3) between SCUBA's commissioning in 1997 
and retirement in 2005 are considered for analysis 0. SCUBA uses 
two jiggle mapping modes: 16 point, which fully samples only the 
850 /xm array; and 64 point, which fully samples both the long and 
short wavelength arrays. Either sampling mode can be used here, 
although because the 450 /^m array data is used as an atmospheric 
monitor, data sets where the data from the short wavelength array 
were not recorded cannot be used (this only occurred in approxi- 
mately 1 per cent of the data). Those data sets which were flagged 
as 'aborted' during the observation may have major defects (incor- 
rect telescope pointing, focusing errors, etc.) and are checked and 
rejected if necessary at this stage. A summary of those data sets 
which pass these initial criteria are given in Appendix A. 



2.2 New C10152. 7-1357 observations 

The work of lZemcov et al.l l l2003h showed that heterogeneous chop- 
ping patterns make SZ effect measurements with SCUBA very dif- 
ficult. To test our ideas about optimizing SZ effect measurement 



'http : / / cadcwww . hia . nrc . ca/ jcmt /I 
^ The only clusters which have been excluded are those studied as part of 
the Red Cluster Survey l ,Gladders & Yee..2Q05.) : these cuiTently have little 
supporting SZ decrement and X-ray data, making it difficult to construct a 
model for their expected SZ effect brightness distribution. 
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with a spatially chopped instrument like SCUBA, new observa- 
tions of a high redshift galaxy cluster were performed using a 
chopping strategy designed to optimize the signal from the clus- 
ter's SZ effect increment. These new JCMT data were obtained 
on Sept. 28, 29 and Oct. 17, 2003. The target of these observa- 
tions was CI 0152.7—1357, the most distant cluster in the ROSAT 
galaxy cluster catalogue. This particular cluster was chosen for new 
observations because it lies at high redshift and, due to its X-ray 
brightness, is well studied at other wavelengths. CI 0152— 1357 
was originall y discovered in the Wide Angle ROSAT Pointed Sur- 
vey (WARPS; [Scharf et aljl997h. and independ ents in the ROSAT 
Deep Cluster Survev (RDCSi lRosati et aljl998l) and Serendipitous 
High-Redshift A r chival Cluster (SHARC) samples jRomer et"al] 
I2OO0I) . IJov et alj l l200lh have previous ly detected the SZ e ffect 
decrement in this cluster. More recentlv. lMaughan et al. I (iool) ob- 
served this cluster with Chandra, and found that it is resolvable into 
northern and southern 'clumps' of X-ray emission. The southern 
clump's X-ray emission is highly peaked, which makes unambigu- 
ous determination of the SZ effect in this cluster difficult. There- 
fore, our observations focused on the northern cl ump, centred at 
aj20o o = l''52'"44.=2, 5j2ooo = -13°57'16" l lMaughan et all 
l2003h. The data were taken in a similar manner to those discussed 
in IZemcov et alj j2003h . A 16 point jiggle pattern with multiple 
chop throws and orientations was used to sample the spatial region 
around this cluster. As chopping onto the southern clump would be 
highly undesirable, the chops were carefully selected to avoid this 
region. The average opt ical depth as measur ed by the JCMT water 
vapour monitor (WVM. Iwiedner et alj20ofl) during these observa- 
tions was rgso = 0.225. In all, 13.4 ks of data were taken for this 
field. 



2.3 Preliminary analysis & calibration 

The data are first double-differenced and flat-fielded using the 
standard SURF analysis packagfl Double-differencing involves 
subtracting the data between the two nod positio ns to remove both 
long term drifts in the data jHoUand et alj|l99^ an d a systematic 
effect which may be due to differential illumination jZemcov et alj 
I2OO5I) . Fortunately, none of the data discus sed here s uffer from the 
16 sample correlated noise effect noted bylBorvs et a l. (2004 ) and 
subseq uently discussed by IZemcov et all 1 20051) and IWebb et alj 
( l2005l) . so we can proceed without correction. The flat-fielding pro- 
cedure multiplies the data from each bolometer by a known amount 
to account for variation between detectors; the multiplicative val- 
ues are very stable in SCUBA data (T. Jenness, private communi- 
cation). 

Next, the SURF package is used to correct the data for atmo- 
spheric extinction. This involves multiplying each bolometer's time 
series by e^^''""^, where A is the airmass of the pointing, rgso is 
found from either the CSO r-meter data provided in the JCMT 
archive, using the relationships between T220 ghz and ghz pro- 
vided in [Archibald et alj ( l2002h . or archival skydips should the r- 
meter data be unavailable. The data from each bolometer are then 
de-spiked by removing points greater than 5a away from the mean 
of the time stream. 

Because these data sets sometimes involve observations of the 
same field performed under very different conditions with differ- 
ent observational strategies, in some cases separated by as much 



as 5 years, calibration is a concern. The data for each observation 
are calibrated at this stage by multiplying each time series by a flux 
conversion factor (FCF). For every night in our data set, we retrieve 
all the calibration observations of Mars, Uranus, CRL 618, HL Tau 
or CRL 268 8, which a re the best non-variable calibrators available 
to SCUBA djenness et al. 2002 ). An 850 ^m FCF is derived from 
ea ch of these observati ons, using the peak flux method discussed 
in Ijenness et alj ( |2002|). Planet brightnesses are derived using the 
FLUXES package JPrivett. Jenness & Matthewsl 19981) , whil e the 
secondary calibrators' fluxes are given in Jenness et alj ( I2OO2I) . 

The FCF chosen for each data set is generally that derived 
from the calibration observation closest in time to the observa- 
tion of interest. However, because focusing the telescope tends to 
change the calibration somewhat (T. Jenness, private communica- 
tion), a calibration observation performed right after a focusing is 
not applied to data taken before that focusing. This helps to ensure 
that data and calibration were both taken with the same optical con- 
figuration. In cases where no calibr ations are available, the fiducial 
FCFs given in ljenness et al. H2oo3) are applied to the data. The set 
of calibr ations used here hav e statistics very similar to those pre- 
sented in ljenness et al. I( l2002t) . 



^ For details on the implementation of these procedures, see 
|http : / /www .starlink.rl.ac. uk/star /docs/sun2 16 . htx/sun2 



2.4 Residual atmospheric emission removal 

Bolometers in the SCUBA array each have their own noise prop- 
erties, and these must be checked from night to night. As is com- 
monly done, we remove the noisiest bolometers in the array from 
analysis. For each data file, individual bolometers with time-stream 
variances larger than 77 times the average bolometer variance at a 
given wavelength are removed from the analysis pipeline; in this 
analysis, 77 = 1.5. While this value of 77 may appear to be surpris- 
ingly low, because the very noisiest bolometers dominate the vari- 
ance 7; ~ 1 does not imply we are removing half of the bolometers. 
Typically only between 2 and 7 of 37 bolometers for a given data 
subset are removed for this 77. Although we could perform an it- 
erative procedure (with a larger value of 77) to remove bolometers 
based on a recalculated average variance, we find that this is unnec- 
essary. 

Thermal emission from both the atmosphere and telescope's 
surroundings contribute a large spurious signal at sub-mm wave- 
lengths. The double-differencing scheme SCUBA employs reduces 
much of this signal, but atmospheric effects still persist, mainly as 
a common-mode noise across the array. It is standard practice to 
subtract the array average at eac h time step to remove these resid- 
ual signals dHolland et alj|l999l) . Unfortunately, the SZ profile is 
appreciable compared to the array size. This means that removing 
the 850 fim array average at each time step would remove a non- 
negligible amount of the SZ signal in these data sets. Fortunately, 
the atmospheric no ise is strongly corre l ated between the arrays at 
both wavelengths ([Holland et alj 1 1 999l iBorvs. Chapman. & Scotj 
Ijenness et alj l2002h . Therefore, we can use the 450 /im time 
stream to subtract the common-mode atmospheric contribution 
from the 850 fim time stream; this is accomplished as follows. 
Denote the time series at wavelength A as Rt{t), where A G 
{450 ^im, 850 fxm}, b is an index running over bolometers and t is 
time. The 450 ^m and 850 /im array averages, denoted are 
calculated at each time step. Because double-differenced SCUBA 
data retains some arbitrary instrumental offset in its time series at 
iStham a avelength, we also calculate the DC offset of both time se- 
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ries, R^, via 

w=^Y.R\t), (1) 

where T is the length of time in the data subset under considera- 
tion. The data series M*^°{t) = {R'^^°{t)) - R^ is then formed, 
and is fit to {R^^'^{t)) to obtain a multiplicative scaling factor, a. 
An atmosphere subtracted 850 fim time series, Sf ^'^ (t), is then cal- 
culated via: 

Sf " ) = iif ) - a (ii*'^'' ) ) - R^, (2) 

as disc ussed in lBorvs. Chapman. & Scotl ( 1 19991) and lZemcov et al.l 
( I2OO3I) . This Sl^°{t) is the final, cleaned time series used in the 
analysis. The value of a does not vary by more than 10 per cent 
over a nig ht, but does vary considerab l y on t ime sc ales of months. 

Both iBorvs. Chapman. & Scot^ Il999l) and [Archibald et"ai] 

( I2OO2I) discuss the correlation between the bolometer-averaged 
time streams of SCUBA's 450 fim and 850 nm arrays in some de- 
tail. However, the effectiveness of using the 450 fim atmosphere 
removal method for large data sets has not been considered in pre- 
vious work. As such an investigation is critical to understanding the 
errors in our SZ estimates, we perform a new analysis here. 

The variance of the time stream of each filqj is calculated af- 
ter double-differencing, extinction correction, etc., but before any 
array average has been removed; this is termed the 'raw' variance, 
""raw Atmospheric removal using the array averages at each time 
step is then implemented using either the 850 /im amy average 
(i.e. the standard StJRF) method, or the fit to the 450 /im array 
average method (Eq.[2l(. The variances of these 'atmosphere sub- 
tracted' time streams (denoted by o^^r^Q and a'L^rM' respectively) 
are calculated, and the variance ratios g^-gso ~ o"raw/o"-850 ^nd 
g?.450 = o"?aw/f-450 founcjf]. Fig. [T] shows for a random 
sample of data for the two atmospheric subtraction techniques; it is 
worth noting a few characteristics of these plots. 

The absence of points below = 1 suggests that subtract- 
ing the residual atmospheric via either method never increases 
the noise in the time stream, and therefore should always be im- 
plemented. This subtraction removes residual 1// noise in each 
bolometer's time series left after double-differencing. Subtracting 
the scaled 450 /im average is never as effective as subtracting the 
850 lira array average, however, it is not a drastically worse tech- 
nique, as the average ratio g-85o/9-450 is about 0.9. Interestingly, 
there seems to be little correlation between the weather (as mea- 
sured by '^850 ^ the airmass of the observation) and the improve- 
ment in variance. 

Because it is critical that the analysis method preserve the 
large scale structure in the data, the 450 /im subtraction method is 
applied in our SZ effect pipeline. However, as using the 850 /im 
data is more efficient for removing atmospheric noise, this ap- 
proach can be used when searching for point sources in the maps. 
Therefore a pipeline has been created which uses the 850 ^.m aver- 
age subtraction method to reduce data for finding point sources in 
the maps (as we discuss in the next sub-section), while the 450 /im 
atmospheric subtraction method is applied to the data for the SZ fit. 

* Here, one file consists of between 20 minutes and an hour's worth of 
data; the exact duration of these files is unimportant, as the central limit the- 
orem applies for anything more than a few minutes of double-differenced 
SCUBA data (> 100 samples). 

^ It is important to note that is the average variance of the bolometers' 
time series, not the variance of the average time series. 



2.5 Point source removal 

After atmospheric removal, the 850 /im data are binned into a sky 
map with a pixel size of 3 x 3 square arc seconds. Any pixel ob- 
served less than 20 times is removed from fu rther analysis, as th e 
noise in these can be highly non-Gaussian jCoppinet aT]|2005h . 
This procedure generally rejects less than I per cent of the pixels in 
a map, so it is not a very costly cut. In mosaicked fields involving 
multiple pointings, the data are simply co-added. Removing noisy 
bolometers, calibration and performing atmospheric removal be- 
fore this binning ensures that the data in such mosaicked fields are 
treated on an equal basis. A standard deviation map is also made 
for each field; this provides an estimate of the error in each pixel. 
These two maps are convolved with a Gaussian point spread func- 
tion (PSF), resulting in a flux map and an error map for each field. 
The flux map is divided by the error map to produce a signal to 
noise ratio (S/N) map, which are shown in Fig.|2] The source fluxes 
and their significance in these maps are the same as those one would 
find by fitting the JCMT's PSF to unconvolved maps. The maxima 
in the flux or S/N maps are the locations of possible point sources. 

The method used to identify possible point sources in these 
fields consists of selecting those points in the S/N map above a 
given threshold. For short integration time maps, which we de- 
fine to have total integration times tint < 15 ks as listed in Ap- 
pendix A, this threshold is 4(7, while for long integration time maps 
tint ^ 15 ks, we have chosen the threshold to be Scr. These limits 
are chosen because the short integration time fields have a much 
higher incidence of spurious sources at the 3(t level than would be 
expected based on the quoted noise levels. A Gaussian PSF of the 
same flux as the removed point source is then subtracted from the 
unconvolved map, and the for the difference is calculated. This 
statistic provides a measure of how similar the point source is to the 
telescope's fiducial beam; the area around a peak in the map poorly 
fit by the telescope's PSF is likely due to random fluctuations, addi- 
tional fainter sources, or systematic noise effects, rather than a real 
source. 

In some cluster fields significant point source detections oc- 
cur where the SZ signal is expected to be strongest, within about 
30 arcsec of the measured X-ray centre (these sources are noted in 
Appendix B). As naively subtracting flux from this region would 
erroneously decrease the measured SZ effect increment, we do not 
subtract these central sources from the time series. This can be jus- 
tified by considering the characteristics of the sources of sub-mm 
emission in these cluster fields. 

Gravitational lensing is known to magnify background sub- 
mm sources in many cluster fields. Lensing models de-magnify 
sources in the interior of the Einstein ring, even when relatively 
complex mass distributions are considered. Sub-mm sources be- 
hind the cluster would either be lensed to the Einstein ring and 
beyond, or be much dimmer than their unlensed brightness, thus 
leaving the central area of a cluster free from lensed emission. This 
type of source would only be problematic in low mass clusters, 
where the Einstein ring is close to the centre of the cluster. Fortu- 
nately, the SCUBA catalogue almost uniformly contains massive, 
rich clusters, so we expect little contamination from these sources. 

The possibility of dust emission from the central galaxy is 
more problematic. As no measurements of pure dust emission in 
the far-infrared (FIR) in mid- to high-redshift cD galaxies exist, 
it is diffi cult to estimate the flu x contribution expected from these 
sources. IVlahakis et al.l | |2005|) show that early-type galaxies of a 
given mass are dimmer than comparable late-type galaxies, and 
since even canonical Ultraluminous Infrared Galaxies (ULIRGs) 
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Figure 1. These plots compare the atmospheric removal methods. The left panel shows g?_g5o ^ function of weather, as traced by ^ X rgso (where A is the 
airmass of the observation) for a random sample of 250 data files in our set. The right panel shows the quantity a'?_45o versus A X Tgso for the same random 
sample. The histograms to the right of either of these plots show the distribution of both q^. The ratios are always greater than 1, implying that some form 
of atmospheric removal should always be instituted. However, improvement from the 450 fim atmospheric removal method is systematically lower than for 
the standard 850 fim method. The average of the ratio q^g^fy/q'LigQ is about 0.9, meaning that the difference between the improvements for the two methods 
is not substantial. The overall improvement in variance after atmospheric removal is typically about a factor of 2; we find that generally, removing an array 
average does not improve an individual bolometer's noise performance by a large factor However, some sort of array average subtraction is critical for the 
noise characteristic of the time series, as it removes residual 1// noise left after double-differencing. 



have fluxes less than ImJy ai z = 0.5, a priori we expect little 
contamination from typical central cluster galaxies. 

Should they in fact be sub-mm bright, there are two meth- 
ods which could be used to identify central galaxies in these clus- 
ters. The first relies on the fact that, if present, such sources would 
be very bright at 450 /im. It is therefore necessary to check the 
450 nm maps for sources coincident at both SCUBA frequencies; 
these would indicate that the central flux is indeed due to bright 
dust emission in the cluster's central galaxy. For example, con- 
sider the case of a 2 — 0.5 850 /im central source with flux 5 mjy, 
similar to the situation in CI 0016-1-16. A conservative estimate for 
the 450 fj.m flux of such a source based on dust emission alone is 
~ 40 mJy, which would be detected in the 450 ^m maps. Further- 
more, in the class of clusters with detected central 850 /im emis- 
sion, C10016-f 16 has the second highest redshift, meaning that 
the 450 fim flux of central sources in the other cluster fields would 
be much stronger stilQ The second method of identifying central 
cluster sources involves measuring the shape of the source in the 
map. A number of the central sources in these clusters appear to 
exhibit extended structure poorly fitting the telescope's PSF Such 
sources (as in Abell 370, Abell 383, Abell 520, Abell 1689, and 
RXJ1347— 1145) are unlikely to be galactic, as these should be un- 
resolved. 

This being said, a counter-example to the generalization that 
moderate redshift cD galaxies are sub-mm dim is the central 
galaxy detected in the z = 0.25 cluster Abell 1835, which multi- 
wavelength studi es have shown to be an FIR bright cD galaxy 
( lEdge et al.|[T99^ . The image of this source in the SCUBA map 

^ The highest redshift cluster with a centrally-detected source, 
C11604-I-4304, is a poor choice for SZ effect measurements using 
the methods described in this work, as discussed in Section [3.2.1| 



shows it to be a point source well matched by the JCMT's PSF 
rather than extended emission, so even in the absence of measure- 
ments at other wavelengths it would be considered suspicious. This 
source is the only sub-mm bright cD galaxy identified in the liter- 
ature and detected in our sample, and we argue that it is an excep- 
tional example. In any case, this field is flagged as contaminated 
for the purposes of SZ effect measurement with these data. 

It is certain that central cluster galaxies emit at 850 nm at some 
level; if the average brightness of these galaxies is ~ 1 mJy, which 
is much brighter than we would naively predict, they would be a 
significant contaminant for measurements of the SZ effect (this is- 
sue is discussed further in Sections 13.31 and \4\. The sub-mm flux 
distribution of cluster sources is certainly an open topic, and one 
which merits further study. 

Table [5] in Appendix B lists the point source candidates, po- 
sitions, fluxes and their errors, and the statistics for the binned 
maps. As none of the unresolved sources exhibit a poor spatial 
statistic, they all meet our criteria for being real. For the purpose 
of measuring the SZ effect signature, all detected point sources 
must be removed from the data. We do this by subtracting all of 
the sources in Table |5] directly from the time series by sampling 
their idealized shape (i.e. the PSF shape with the same amplitude 
as the source's brightness) and subtracting the resulting time se- 
ries from the real data. These cleaned data are used in subsequent 
stages of the analysis. Appendix C provides notes about each of 
these clusters, including references to previous SCUBA measure- 
ments of these fields. 

In order to provide a check for the point source identification, 
the 450 fim data are also reduced using the SURF pipeline. The esti- 
mated 450 fim fluxes at the positions of the 850 ^m sources are also 
listed in Appendix B. Only one 850 ^m source is robustly detected 
at 450 fim, RXJ 1347-1145-1. 
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Figure 2. 850 signal to noise ratio maps of clusters in the SCUBA archive. The signal to noise ratio is used because the noise chai'acteristics of these 
maps can be quite heterogeneous. Note that these images are useful for assessing the quality of the data and the contamination by bright sources, but do not 
show the entire SZ effect (since a large part of the signal comes from differences relative to the off-beams). The cross symbol shows the centre of the cluster 
as defined by X-ray data. For the long integration time fields (tint 15 ks), the dashed contours are at {—7, —5, —4, —3} X a and the solid contours are 
at {3, 4, 5, 7, 9, 11, 13} x a, while for the short integration time fields (4int < 15 ks), the dashed contours are at {—7, —5, —4} X a and the solid contours 
are at {4, 5, 7, 9, 11, 13} x a. 'Holes' in the maps occur due to noisy bolometer subtraction, and generally occur only in fields with short integration times. 
Negative beams of sources can be seen in some of the maps where 'in-field' chopping was used. Abell 370 exhibits a bright source near the southern edge of 
the field, Abell 478 exhibits a bright quasar near the middle of the field, and Abell 496 shows emission associated with its central galaxy. 
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Figure 2 - continued Abell 780 contains the bright source Hydra A, which is visible in this image, which, due to short integration time is sparsely sampled. 
Abell 1689 contains an extended source near the X-ray centre; we associate this with the SZ effect (see Section|4). Abell 1835 also shows a significant detection 
at its centre, although this is associated with abnormally bright dust emission in the central cluster galaxy. Abell 2163 has two pointings, neither of which are 
very deep, and both of which show patho logies in the data; the cluster is removed from further analysis because of this. Abell 2218 contains a multiply imaged 
high redshift source iKneib et alj|2004hl) ; since this would interfere with SZ effect extraction, this cluster is also removed from further analysis (because the 
central source in this cluster is very significant, the contours in the Abell 2218 map are drawn at (—5, —4, —3, 3, 6, 9, 12, 15, 18, 21) X cr). 
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Figure 2 - continued Abell 2597 contains a bright poi nt source associated with AGN activity in its central galaxy. CI 0016+16 has a central enhancement 
which we ascribe to the SZ effect lZemcov et alj )2003h . CI 0055— 2754 has a very short integration time and so is sparsely sampled. CI 1604+4304, a high 
redshift cluster, appears to contains a large number of sources. 
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Figure 2 - continued CI 2129+0005 is anothe r field with an ex tremely short integration time. MS 0451— 0305 contains another resolved gravitationally 
lensed source (or set of sources) at high redshift jBorvs et al.l2004 . Because the map of RXJ 1347—1145 covers such a large area, it appears twice; the lowest 
left map is the whole field, and the lowest middle map is a detail of the centre of the cluster. The central area of RXJ 1347—1145 contains a source bright at 
both 850 and 450 ^m, and, potentially, a resolved image of the SZ effect in the centre of the cluster. 



3 SZ INCREMENT FITTING 



3.1 Field cuts 



Although all of these clusters should exhibit the SZ effect, the sig- 
nal from some may be swamped by known sub-mm bright central 
galaxies, complex emission due to gravitational lensing, and in- 
strumental pathologies. Therefore we cull the list to exclude those 
fields where SZ emission at the ~ 1 mJy beam^^ level cannot pos- 
sibly be disentangled from other sources of emission based on our 
present knowledge of the cluster. 

We remove those clusters lying at low redshift (2 < 0.1) with 
bright AGN near their centre which would swamp the SZ effect. 
The clusters affected by this cut are Abell 478, Abell 496, Abell 
2597, Abell 780, and CI 2129+0005. In general, these clusters 



were only observed for a short time, and have poor spatial coverage. 
Similarly, Abell 1835 is removed from the sample as it contains an 
abnormally bright cD galaxy coincident with the X-ray centre, as 
discussed in the last section. 

A further cut re moves two clusters w ith strong lensing sig- 
natures: Abell 22 18 jKneib et all l2004bl) and MS 0451-0305 
( iBorvs et al]|2004h . In these systems, it is too difficult to disentan- 
gle the complex lensed structure from the SZ effect. Abell 1835 is 
also cut from t he sample, as it is known to harbour a sub-mm bright 
central source tedge etal.ll999l) . 

The two clusters (CI 0055- 2754 and CI 2129+0005) with 
under 30 minutes of integration time are removed from further con- 
sideration; the noise level achievable in this time under even the 
best conditions is too high to see a 1 mJy beam^^ SZ effect. 
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Furthermore, neither X-ray data nor velocity dispersions are cur- 
rently available, so we cannot apply our SZ fitting method to these 
clusters (as discussed in Sectior[3]2j. 

The noise structure in the map of Abell 2163 is atypical, 
demonstrating elongated patches and deep negative holes. Despite 
our extensive background with SCUBA data, we are not able to 
explain this noise pathology, and therefore flag this cluster map as 
suspect and remove it from our sample. 

After these cuts are applied, our initial SZ effect sample con- 
sists of 33 galaxy clusters. Using the version of the SCUBA time 
series and maps with point-sources removed, we now fit these data 
for the SZ increment. 



3.2 Sunyaev-Zel'dovich effect modeling & fit 

The isothermal /3 model is commonly used to parameterize the 
spatial distr i bution of the SZ effect, as discussed in, for example, 
lBirkinshawl ( ll999h . If is the angular distance from the centre of 
the cluster, I{6) is the SZ effect intensity at 9, Iq is the central SZ 
effect intensity, and /? and 6c are values parameterizing the cluster, 
then the isothermal (3 profile is given by 



1(0) = -fo 1 + ^ 



2\ (l-3/3)/2 



(3) 



Given the low signal to noise ratio of our SCUBA data, we cannot 
hope to determine all of these parameters from our data alone. For- 
tunately there is a rich body of ancillary data that can be used to 
assist in the fits. 



3.2.1 Independent measures of 9c and (3 

28 of 33 clusters in our sample have published X-ray and/or SZ 
decrement measurements from which 9c and /3 have been deter- 
mined. These parameters and references are listed in Table[T] 

Of the remaining 5 clusters, two (Abell 222 and Abell 851) 
have archival X-ray datfl which we use to fit the center coordi- 
nates, dc, and (3 ourselves. To ensure our fitting procedure was ro- 
bust, we downloaded and fit the model to archival data for clusters 
with already published 9c and (3 values, and found no significant 
discrepancies. 

For the three remaining clusters for which we could find 
no archival X-ray data (CI 0303+1706, CI 1604+4304 and 
CI 1604+4321), we assume the clusters are isothermal and have 
a standard King profile (i.e. /3 = 2/3). An estimate of the core ra- 
dius rc (which is related to 9c by the angular diameter distance) can 
be obtained via t he line of sight velocity dispersion of the cluster 
Gr l lSaraziiJl98^ : 

,2 



(4) 



Here a value for po must be adopted, since no direct measure- 
ments of the central density of these clusters exist. If we esti- 
mate rieo = 10^^ cm^'' and assume that globally the intra-cluster 
plasma is electrically neutral, then 



po — Ueoime + nip), 



(5) 



where nie and irip are the electron and proton masses, respectively. 
This calculation yields the core radius estimates listed in Table [T] 



^ Reduced maps from NASA's High Energy Astrophysics Science Archive 
Research Centre were obtained. 



A position must also be adopted for the X-ray and SZ peak. 
For CI 1604+430 4 and CI 1640+4321, ^ye ass ume that the optical 
centre as given in lOunn. Hoessel & Oke 1 19861) corr esponds to the 
X-ray peafl For CI 0303+1706. IUeda et alj ( 1200 ll) give an X-ray 
centre which we adopt. The model parameter predictions for these 
3 clusters are listed in Table[T] 



3.2.2 Fitting for the SZ increment 

With these parameters in hand, we can fit for the amplitude of the 
SZ increment in the clusters, for which two different approaches 
have been adopted. 

In the first method, a model map is made for each cluster us- 
ing its isothermal (3 profile parameters. This map is convolved with 
the JCMT PSF and normalized to have a maximum of 1.0 mjy 
beam^^; this ensures that the fit discussed below yields a consis- 
tently calibrated result. This model map is then differenced in the 
same way as in the SCUBA observation, making a model time se- 
ries. The mean of the model for each subset of data is removed to 
match the subtraction of the i?*^'' from Equation[2]in the real data. 
The data are then fit to the model differences to determine a scal- 
ing coefficient which gives the central increment value. Because it 
directly incorporates our differencing scheme, this fit method ob- 
viates the need to track the correlations which would occur in a 
pixelization scheme. 

The second method for determining the SZ effect amplitude 
involves making a map using the S'f (t) and fitting the result to 
the model. In this method, a model map is made and normalized 
in the same way as discussed above. This map is then differenced 
according to the actual SCUBA observation to produce a time se- 
ries. Again, the 7?*=" is subtracted from the model time series. This 
model data is then rebinned into a map in the same way as for the 
real data. The data map is then fit to the rebinned model time series 
to yield a scaling, which is the SZ effect amplitude for each clus- 
ter. This method allows a check on the robustness of time series 
analysis in each cluster. 

The results of both methods are listed in Table|2l and are plot- 
ted against one another in Figure[3] Although there is some scatter 
between the methods, they generally agree, showing that the fitted 
SZ effect amplitude from these data is largely independent of the 
fit. However, the error bars from the map fit method are often ~ 5 % 
larger than from the time series fit, and so we adopt the SZ effect 
amplitudes given by the first method as the best fitting increment 
values. 

As a test of the assumptions about the effects of central point 
sources discussed in Section |23] these fits are also performed with 
the central 18 arcsec masked from the fit. This test allows a study 
of how strongly the fits depend on the data from the very central 
region of the cluster; an indication that point source contamination 
may be present is if the average fit amplitude is reduced in those 
clusters with central sources (the starred entries in Table |5]l. Since 
data are being removed from each set, the error bars should also 
increase. 

To quantify the change in the fit amplitudes, the ratio of the 
un-masked to masked fit amplitude in each cluster is used; this 
should be unity for no change in the fit amplitude. For the entire 
sample listed in Table|2] the mean ratio is 1.003 and the standard 



° While some archival XMM-Newton data exist for the supercluster field 
containing these two clusters, the data ai'e not sufficiently deep to reliably 
extract model parameters. 
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Table 1. Summary of the isothermal /3 model parameters in the SZ sample. The first part of the table lists the parameters for clusters with good model 
parameters and X-ray centres in the literature. The second part of the table lists model pai'ameters and centres derived from archival X-ray data. The third part 
of the table lists the assumed parameters for clusters with no available X-ray data (how these parameters are dervied is discussed in Section[3j2). 
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° as listed in the NASA/IPAC Extragalactic Database, http : / /nedwww . ipac . caltech . edu/ 

' In order, references are given for the central X-ray/SZ effect position, isothermal /3 model 

parameters and X-ray temperature: (1) This work; (2) Allen & Fabian (1998) : (3) Allen QOOff); 

(4) Allen. Schmidt & Fabian (2002) : (5) Bohiinger et al. (2000) ; (6) Cruddace et al. (2002) ; 

(7) David. Forman & Jones (1999) ; (8) Ebeling etal. (1998) ; (9) Ettori & Fabian (1999) ; 

(10) Ettori et al. (2004) ; (1 1) Gioia & Luppino (1994) ; (12) Gioia etal. (1998) ; 

(13) Govoni et al. (2004) ; (14) Grego et al. ( 2000) ; (15) Hicks et al. (2002) ; 

(16) Hughes & Birkinshaw (1998) ; (17) LaRoque et al. (2003) ; (18) LaRoque et^ (2006) ; 

(19) Mason & Myers (20()ff); (20) Maug han et al. (2003); (21) Mohr'e t al. (200(J); 

(22) Ota etal. (1998) ; (23) Ota etal. (2004) ; (24) Ota & Mitsuda (2C)04); (25) Reese et al. (2002) ; 

(26) Reipri ch & Bohringer (2002) ; (27) Rizza etal. (1998) ; (28) Rosati etal. (1999) ; 

(29) Sand ers. Fabian & Allen (2000) ; (30) Soucail et al. (2000) ; (31) Vikhlinin et al. (2002) ; 

and (32) Xue & Wu (2002) . 

This assumes the standard King profile, /3 = 2/3; see text for details. 

References for the coordinates and velocity dispersions for these clusters are given in: 
(l) Uedaet al. (2001) ; (2) Gunn, Hoessel & Oke (1986) ; and (3) Postman et al. (2001) . 
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Table 2. Summary of fits for the full sample of SZ effect clusters. 



Cluster Fit to time series differences Fit to maps 

(mjybeam^-^) (mJybeam^^) 



Cluster parameters available in the literature. 
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0.68 


2855763 / 2856074 


-0.64 


0.82 


7457 / 7235 


PKS 0645-19 


-1.45 


2.63 


153568/ 153557 


-2.30 


2.97 


2680 / 2778 


RXJ 1347-1145 


8.80 


0.84 


876116/876350 


9.06 


2.10 


5212/5392 


Zwicky 3146 


1.47 


1.01 


208741 /208776 


1.81 


1.16 


2865 / 2807 



Cluster parameters derived from archival X-ray data (columns as above). 



Abell 222 -24.6 8.0 158325/ 158329 -25.8 8.7 2462/2304 

Abell 851 -10.1 2.5 1328111 / 1327919 -7.24 3.75 4088/4024 



Cluster parameters derived from velocity dispersions (columns as above). 

C10303+1706 5.66 6.29 220745/220752 2.79 6.15 2360/2482 
CI 1604+4304 4.53 0.91 714956/714771 4.34 0.95 3026/2795 
CI 1604+4321 2.51 1.31 242345/242337 3.09 1.36 2466/2453 



deviation of thie ratios is 0.22; tliis means that, on average, the fit is 
unchanged, although there is some scatter in the ratio on a cluster 
by cluster basis. The average ratio for clusters with detected cen- 
tral point sources is 0.995 ± 0.240; this is also consistent with no 
change, on average. In addition to these amplitude changes, the av- 
erage enw bar for the fits increases by 7 per cent, as expected. 



Based on this analysis, we conclude that there is no significant 
difference between the fits for clusters with and without detected 
central point sources. Furthermore, this test suggests that the meth- 
ods we employ are mostly sensitive to the shape of the SZ effect in 
the data, rather than the brightness of the centre-most region alone. 
In the absence of results supporting the proposal that these fits are 
strongly contaminated by central point sources, the fit results which 
include the central regions for each cluster are quoted hereafter. 



3.3 Sunyaev-Zel'dovich effect limits 

Because the typical SZ effect in these clusters is expected to be 
~ 1 mjy beam^ ^ , cluster fields with short integration times will 
have noise levels much greater than the expected signal, and merely 
increase the noise in the overall measurement. Because they pro- 
vide relatively uninteresting estimates of A Jo, the cluster fields 
with low integration times (tint < 15 ks as discussed in Section 
12.5b are excluded from the results presented in this section. Also, 
clusters whose model parameters are derived from the cluster's ve- 
locity dispersions or poor X-ray data are suspect; although these 
may be a reasonable estimate of the model, the potential error in 
the model dominates the statistical error in the data. These clusters 
are also excluded from this section. These criteria leave a final SZ 
sample of 17 clusters, listed in Table[3] 

At 450 /im, the SZ distortion is a small but non-zero fraction 
of its peak. Because the 450 /im array averages, which contain the 
differenced SZ effect flux, have been removed from the data as part 
of the atmospheric removal method, the fits to the isothermal 13 
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Figure 3. Comparison of fitting the SZ effect model to the time series data 
(shown on the abscissa axis) and the binned map data (shown on the ordi- 
nate axis). The points show the amplitudes from the fits, with their associ- 
ated 1 cr error bars. The points with dark centres highlight the 17 final SZ 
effect clusters defined in Section 13.31 The dark dashed line shows the line 
X = y, while the light dotted line shows the result of a linear fit to all 33 
points, and the dai'k dotted line shows the result of a lineal' fit to the 17 final 
SZ cluster points. The fits have slopes very close to 1, and so are difficult 
to differentiate on this plot. This shows that the two methods yield very 
similar results in most of the clusters. The point near {—6, 6} is for Abell 
2219; it is unclear why it lies so far from the expected relation, although the 
appUcation of a poor model is the most likely cause. 



model mus t be corrected for the flux removed by this procedure (see 
IZemcov et al. 2003). The relation between the quantity measured 
by our isothermal [3 fits, A/J*, and the true peak SZ flux at 850 /im, 
is given by 



(6) 



where AI^^^ is the flux of the SZ effect removed by the subtraction 
of the scaled 450 ^tm array averages. To determin e A/ave ^ we adopt 
the procedure discussed in lZemcov et alj | |2003|) . For each cluster, 
the SZ shape as defined by the isothermal /3 model is sampled in the 
same way as in the experiment, and a linear coefficient ^ relating 
the ratio of the input SZ amplitude to the mean of the sampled data 
stream is found. This gives Alt^° = C A/(J^", where ^ includes the 
gain difference between the 450 /im and 850 /im channels. The ^ 
coefficients are given in Table[3] 

For each cluster, the best fit AJg^" and error bar ctsso are 
found using a likelihood function and determining its maximum 
and associated 68 % confidence limits. The likelihood function is 
created by making a set of SZ effect models with different Comp- 
ton yo parameters. At 850 /im, non-thermal corrections to the ther- 
mal SZ effect can change the best fitting amplitudes by as much 
as 15 per cent, given the temperatures of the clusters in the sam- 
ple (although typically the corrections are much smaller, around a 
few per cent for most clu sters). To account for t his, the corrections 
to fifth order discussed in lltoh, Kohvama & Noza wa ( 1998|) are ap- 
plied to these fits, using the electron temperatures listed in Table 
[T] These temperatures are derived from averages over X-ray data 



Table 3. SZ effect parameters derived from ai'chival SCUBA data. Listed 
are the cluster name, correction factor for the 450 /im atmosphere subtrac- 
tion 5, D = A/*f /A/q^*^, and relativistic effect corrected amphtudes in 
both MJy sr~ ^ and (unitless) yo . 



Cluster 




V 


A/8'^o 


yo 






xlO^ 


MJy sr~i 


xlO* 


A Kail TflQ 


U.U44 


1 1 
l.j 


n 99+0.80 

'J-^^-o.so 


1 00+5.02 

i-3»_5 02 


Abell J /U 


0.164 


4.4 


r, r,r, + 0.09 

^•^^-0.09 


1 no + 0.56 

J^-^^-o.se 


Abell J6J 


0.073 


2. 1 




9 20 


A u.^11 con 
Abell 52tJ 


0.255 


7.6 


„ _ +0.17 

0.77_o 17 


A Qn+1.06 

4-82_i 06 


Abell JoD 


U.J J J 


1 c r\ 
1 J.U 


n fi9+0.15 


oc+0.94 
d-»6_Q 92 


Abell 773 


0.006 


0.2 


n on+0.59 

-0-80-0.59 


r- in+3.74 
0-l^_3.74 


Abell 963 


0.043 


1.1 




-0.521^7! 


Abell 1689 


0.277 


10.4 


0.87i«:i^ 




Abell 2219 


0.084 


3.1 




7 n9+3-''6 


Abell 2390 


0.149 


4.8 


o.7il°:ll 


4.00_Q g4 


CI 0016+16 


0.00 


0.0 


0.58l«:i^ 




CI 0024+1652 


0.117 


2.8 


n 04+0.26 




CI 0848+4453 


0.070 


1.5 


-0 ns+'^-"9 

o-i5lS:l^ 


-o.2oi|;:g 


CI 1455+2232 


0.339 


8.2 




MS 0440+0204 


0.142 


3.6 


o-o3l°:i^ 




MS 1054-03 


0.018 


0.5 


17+0-11 




RXJ 1347-1145 


0.135 


4.5 


1 50+0-16 


Q 09+0.96 



for each cluster measured on angular scales similar to those mea- 
sured with SCUBA; this reduces or eliminates the effect of biases 
associated with cooling flows in the central regions of the clusters. 



The difference AJn 



-CA/o- 



is found from this model. The dif- 

D 



ference which most closely matches the observed AJq yields the 
best fitting yo- An integral under the likelihood function is used 
to determine the 68 % confidence limits, which are quoted as 1 cr 
errors in the two values. The best fitting A/q^" is converted from 
mjy beam^^ to MJy sr^^ by multiplying by the integral over the 
Gaussian JCMT PSF with FWHM 14.7 arcsec (equivalent to a con- 
version factor ImJy beam"^ = 5.739 MJy sr"^). The AJo''^°, 
2/0, and their associated errors for each of the final SZ clusters 
are quoted in Table [3] Also included in this table is the parame- 
ter V = ll^'V/o^", which is the ratio of the flux removed due to 
the use of the 450 ^m channel as an atmospheric monitor to the to- 
tal 850 iJ,m SZ flux. This ratio is typically less than 5 per cent, but 
ranges between and 15 per cent in this sample of clusters. 

A common technique used in to draw a significant measure- 
ment from a set of noisy astronomical data is to 'stack' the noisy 
images on top of one another centred on a nominal source position. 
Using the SCUBA maps from the 17 deep integration fields, we 
have performed this procedure with these data. The (non-central) 
sources in these maps are subtracted out before the co-addition 
to give a clean signal from the SZ effect alon^fl These individual 
maps are weighted using the expected central Compton parameter 
in each cluster, which is give n by yp oc \/ilM + (1 + z)'^Q,aTc^'^ 
(as discussed in, for example. iBenson et alj2004|)- He re we use the 
WMAP concordance cosmology jSpergel et al.l2006l) and the clus- 
ter parameters listed in Table [T] In the cluster co-addition, each 



^ If this subtraction is not performed, the co-added map shows both the 
central extended emission am! a diffuse ring of emission with radius slightly 
larger than half an arc minute. This is good evidence that, statistically, dim 
gravitationally tensed sources tend to be re-imaged at or beyond the average 
Einstein ring in these clusters. 
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pixel is given a weight based on thie variance map, and a total vari- 
ance map is made using the same cosmological weight as above. 
The total number of hits in each pixel is also calculated, and this 
can be used to apodize the co-added maps if desired. The resulting 
signal map is shown in Figure|4] 

The co-addition is quite successful at drawing out extended 
central structure in these maps. The peak increment value (even 
without adjusting for the SCUBA chop) is detected at 9a" based 
on the set of 17 deep maps. Comparison of this co-addition with 
a model map is desirable; this is achieved as follows. Using the 
model parameters given above, an isothermal /3 model is created 
for each cluster. This input model map is then convolved with the 
JCMT's PSF and differenced according to the real data. Time se- 
ries averages are subtracted as in the SZ effect fitting pipeline, and a 
model SCUBA map is made using these model data for each clus- 
ter. These model chopped maps are then co-added using the real 
noise weights for each pixel in each cluster. The coadded maps for 
both the data and the model can then be averaged in radial bins to 
facilitate a comparison of the average flux as a function of radius, 
as is shown in Figure[5] 

This comparison shows that the model and data stacked maps 
agree very well from about 0.3 < r < 3 arcmin, where the data 
begin to sparsely sample the sky. However, the first bin (which has 
radius equal to a SCUBA beam) is approximately a factor of two, 
or about 0.75 mjy, higher than the SZ-only model would predict. 
This region is typically not resolved in prior SZ effect measure- 
ments, so these data provide new insight into models for SZ effect 
distributions. The excess emission is difficult to explain in the con- 
text of our current (isothermal /3) model of these clusters, which 
includes only the SZ effect and sub-mm bright sources in the re- 
gion away from a cluster's centre. The extra flux therefore could 
be the signature of unresolved point sources or complex ICM tem- 
perature distributions in these clusters. If this extra flux is wholly 
attributable to the central galaxy in the cluster, a measurement of 
0.75 ± 0.25 mJy is obtained for the average flux of such sources 
over 17 clusters in the range 0.17 < z < 1.2. It is difficult to con- 
sistently apply a more complex ICM model to these data, but some 
fraction of the extra flux may be attributable to varying temperature 
distributions or clumping in the cluster electron gas, as would occur 
in the presence of cooling flows, for example. As the isothermal 13 
model relies on fairly simple physical assumptions (some of which 
are not true in a subset of these clusters), it is quite plausible that it 
fails as an accurate description of the true ICM shape in some frac- 
tion of the sample. However, any other model invoked to describe 
these data would need to provide a factor of 2 greater SZ emission 
inside the central ~ 100 kpc of each of these clusters, or an even 
larger factor in some subset of them. It is most likely that the excess 
emission is due to some combination of a more peaked model, at 
least in some subset of the clusters, and dim galactic sources near 
the clusters' centres. 



4 DISCUSSION 

4.1 Sunyaev-Zel'dovich effect results 

From a rather heterogeneous data set, the methods employed here 
have allowed us to estimate the SZ effect brightness in archival 
SCUBA data for 17 galaxy clusters. For the other 27 clusters with 
archival SCUBA data, we were unable to provide strong constraints 
on the the SZ increment. This is due to a variety of factors, such as 
contamination from lensed sources behind the cluster, AGN in the 



cluster, and low on-source integration times. In general, the com- 
bination of chopping pattern, extended emission and noise means 
that images of the cluster fields within the size of the array do not 
fully represent the SZ signal. However, it appears as though our 
maps of several clusters (particularly Abell 370, Abell 1689, and 
RXJ 1347— 1145) do show significant emission at approximately 
the right brightness and shape to be (chop-filtered) images of the 
SZ effect. 

A number of the clusters observed with SCUBA have also 
been popular SZ effect targets with other instruments (see Ap- 
pendix C for a discussion of this). It is instructive to compare the 
results given here with those of other groups; we have searched the 
literature for other SZ effect measurements for each of these clus- 
ters and Fig.[6]shows the results. 

Comparing to the other SZ effect measurements in these clus- 
ters shows that the SCUBA yo results presented in this work are 
substantially larger than those determined from data taken at SZ 
decrement frequencies in 5 of 8 clusters. Although there appears to 
be a correlation with the results of SuZIE II (Benson et al. 200^, 
which uses data from frequencies on both sides of the null of the SZ 
effect, this is difficult to quantify with the data at hand. The system- 
atic preference for larger yo evident in the SCUBA measurements 
in comparison with measurements at other frequencies may be due 
to some poorly accounted for effect in the analysis or methodology; 
we discuss these here. 

Given these results, a primary candidate for systematically in- 
creasing the measured SZ flux is point source contamination near 
the centre of the cluster fields. Even a relatively dim radio/mm point 
source near the centre of these clusters could serve to increase the 
yo measured from observations taken at SZ increment frequencies, 
and decrease the yo measured from SZ decrement observations. 
On average, such sources would thus increase the disagreement 
between measurements on either side of the SZ effect null. This 
being said, the results presented in Section |3.2.2| suggest that point 
source contamination, at least on average, is not a strong effect in 
these data. Similarly, experiments targeting the SZ decrement are 
typically careful to quanti fy the number and fl uxes of point sources 
in their target fields (e.g. ICoorav et al.lll99§) . On the other hand, 
the excess emission found in the average of 17 cluster maps dis- 
cussed in Section [33] suggests that as much as 0.75 mJy of the cen- 
tral flux cannot be explained by the isothermal 13 model alone. This 
is roughly the same size as the expected SZ effect in these clus- 
ters; we would therefore expect an appreciable bias to affect the 
fits due to central point sources below the detection thresholcPl It 
is not clear whether the excess flux is in fact due to point sources, 
however. 

A second important concern is the model assumed in each 
cluster. The isothermal (3 model, while an accurate description of 
the ICM flux distribution in so me clusters, fails badly in others. For 
example, iBenson et al.l ( l2004h show that the assumed model can 
have a large effect on the yo amplitudes measured in their cluster 
set. To attempt to quantify the dependence of our fitted amplitudes 
on the model, the model parameters are varied and the data are 
reanalyzed. We find that the fitted SZ increment amplitudes are a 
weak function of 6'c: a systematic increase of a factor of 2 in 6c for 
all of the clusters listed in Table |2]leads to an increase in the fitted 
A/o of 15 ± 4 per cent (the quoted error bar reflects the standard 



^" Although not a factor of two, as increasing the value of a single pixel 
in the flux distribution does not have a large effect on the best fitting value; 
this is not true of the of the fit, however. 
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Figure 4. This figure sliows a staclced co-addition of all 17 deep integration cluster fields discussed in Section |33] The weight function used in the stack is the 
expected yo based on the measured Tc in the cluster. The cluster maps are all centred on Aa = 0, A5 = so the imaged SZ effects should all be centred at 
the same point. The image grey scaling is from — 2mJy to 2mJy, and contours are plotted at {—1.0, —0.5} mjy (dotted) and {0.5, 1.0, 1.5} mJy (dashed). 
The statistical noise in this map is \a = 0.2 mJy at the centre of the field, and the map has been apodized using the summed hit map. The box at upper right 
contains an artificial point source with peak brightness equal to that of the SZ effect in this map (1.83 mJy) as a reference for the emission at the centre of the 
co-addition, which in comparison is quite extended. 



deviation of the per cent change for the set). A decrease in 6c by the 
same factor - that is, increasing the 'peakiness' of the clusters by a 
considerable amount - leads to a change in the average j/o of 2 ± 4 
per cent. Alternatively, the fit amplitudes are a stronger function of 
fi: increasing /3 by a factor of 2 leads to a 25±7 per cent decrease in 
the measured A/q. Clearly, a moderate misestimation of the model 
can have an appreciable affect on the ya measurements discussed 
here. 

An issue related to the cluster flux distribution model is that of 
cooling flows: such clusters have more peaked centres than their ap- 
proximately isothermal counterparts. In the set of clusters shown in 
Fig.m Abell 1689, Abell 2390 and RXJ 1347-1145 are classed as 
clusters having cooling cores. Although at least 2 of these 3 clusters 
exhibit SCUBA-determined j/o significantly higher than that mea- 
sured at the decrement, this is not statistically different than for the 
non-cooling flow clusters. More data would be required to reliably 
differentiate between these cluster types. 

Interestingly, Fig. [5] suggests that the centres of our cluster 
sample, on average, emit a factor of 2 more flux at 850 ^m than the 
SZ models used here can account for. This flux may in fact be the 
signature of excess emission from cooling flow core clusters in the 
total cluster sample of 17 fields, or a more generic failure of the 
isothermal (3 model. Equally, it could be the signature of dim point 



sources in these clusters. It is not possible to differentiate between 
these options with the data at hand. 

Finally, two components of the analysis pipeline could have 
a systematic effect on the final amplitudes, the first of these being 
calibration. We have used the SCUBA calibration observations de- 
termined between 2 and 5 times nightly at the JCMT to calibrate the 
observations presented here. Although it is difficult to model how 
the calibrations could be systematically inconect by an apprecia- 
ble factor, there is certainly an approximately 10 per cent variation 
in the calibration runs taken over a single night. For deep cluster 
fields, this effect should average out; however, for shallow fields 
this may not be true. Fortunately, the 8 deep fields discussed in this 
section should not suffer from systematic calibration problems. A 
second component of the analysis pipeline which may systemati- 
cally affect the measured ya is the correction for the 450 /im atmo- 
sphere subtraction. Although the size of the correction is typically 
less than 10 per cent, it may help artificially increase the measured 
yo if the correction algorithm is erroneous. We have attempted to 
be completely consistent with the 450 /im correction, but as it is it- 
self based on the SZ effect model, the correction does increase the 
systematic uncertainty of the measurement. 

As a final compari s on, w e can consider the amplitudes deter- 
mined bv lZemcov etal] | |2003|) for CI 0016+16 and MS 1054-03. 
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Figure 5. Radial averages of the model and data coadded maps for the 17 
deep integration time cluster fields. The diamond points show the averages 
of the stacked data map (shown in Figure |4) in annular bins with width 
15 arcsec (a SCUBA beam). The en'or bars are similarly determined from 
the stacked variance map. The model, shown as a solid line, is determined 
by creating a simulated 'map' for each cluster, sampling it according to the 
telescope's recorded pointing, and running the resulting time series thi'ough 
the data analysis pipeline to create maps which reflect the spatial filtering of 
scuba's chopping pattern and the pipeline itself. These individual model 
maps are then stacked using the same noise weights as are applied to the 
data map stack, and the resulting map is averaged in annuli. To apply the 
proper normalization, the model averages are fit to the data averages before 
plotting. This figure shows that the model and data agree quite well, except 
in the central bin, where the measured flux exceeds the model expectation 
by a factor of ~ 2. 



In that work, the measured yo from the SCUBA data alone were 
(2.2±0.7) X 10"'' and (2.0±1.0) x lO"'*, respectively. Although 
these values are different than those listed in Table |3] the disparity 
is not statistically significant. 

4.2 Point sources 

As part of this survey of archival SCUBA data, we have presented a 
list of point source candidates found in our maps. Broadly, these fall 
into three categories: bright central sources associated with AGN 
activity in the cluster's central galaxy; lensed sources in moderate 
and high redshift clusters; and star forming sources in the clusters 
themselves. Sources belonging to the first category are generally 
very easy to identify, as they correspond to well known galaxies 
at other wavelengths (e.g. Hydra A). The second ca tegory have 
typical ly already been well s tudied, as in Abell 2218 jKneib et al.l 
l2004bl) and MS 0451-0 305 jBorvs et al.ll2004l). and correspond t o 
sources around z ~ 2-3 dSmail et al.lll99HChapman et al.ll2002h . 
Such lensed sources are common in rich cluster fields; most clus- 
ters with long integration times exhibit this type of source. 

Low signal-to-noise sources in a population with steep source 
counts tend to appear somewhat brighter than they actually are. One 
can attempt to correct for this by performing a flux de-boosting 
procedure (e.g. ICoDpin et al.ll2005l) . We have not done this for our 
source list, since the procedure would be complicated for such an 
inhomogeneous data-set, and would also depend on modeling the 
lensing effects. It is therefore important to appreciate that the fluxes 
given in Table |5] will generally be biased high, this being particu- 
larly true for the low significance sources. 



SZ effect is also a strong lensing cluster. This underscores the need 
for experiments with high spatial resolution and multi-wavelength 
coverage, particularly at sub-mm wavelengths. Although the peak 
of the SZ effect and that of maximum lensing have different spatial 
shapes, an instrument like SCUBA - with its single SZ effect obser- 
vation frequency and relatively poor spatial resolution - will have 
trouble disentangling the effects in clusters where the SZ shape is 
broad. 

Nevertheless, the measurements presented here are instructive. 
They show that it is possible to provide useful constraints on the 
SZ increment using mm/sub-mm array receivers and highlight the 
most likely sources of systematic error in any sub-mm SZ effect 
experiment: chopping effects (necessitated by the strongly varying 
atmosphere in the sub-mm regime) and the contamination of point 
sources, particularly lensed ones, which are thought to be corre- 
lated with the existence of a strong SZ cluster. We have shown that 
one must to try to remove the effects of bright point source con- 
tamination in these type of measurements, and that when chopped 
data are used one must be careful to compensate for the effects of 
chopping as much as possible in order not to remove a significant 
part of the signal. 

In the future, the combination of SZ data at a range of wave- 
lengths (from 10s of GHz to 400 GHz or higher), together with X- 
ray, lensing and other optical studies, should allow for constraints 
to be placed on the kinetic SZ effect. This effect will yield the radial 
velocity of the cluster which, in principle, will constrain large scale 
flows in the Universe. Multi-wavelength SZ data will also be use- 
ful for studies of cluster astrophysics, from temperature and mass 
determinations to sub-structure within the cluster gas itself. 

The next generation camera for the JCMT, SCUBA-2 
( [Holland et al.l200llHoiland et al.l2do^ will alleviate many of the 
problems with SZ effect increment measurements highlighted by 
this work. SCUBA-2 will bring CCD-style imaging to sub-mm as- 
tronomy; sampling at 200 Hz, it will allow imaging of the sky much 
faster than atmospheric variation (~ 1 Hz). The large array coupled 
with rapid scanning and map-making techniques should allow for a 
robust estimate of modes on scales many times larger than SCUBA- 
1 has permitted. Since typical increment amplitudes for rich clus- 
ters correspond to ~ 1 mjy per 850 fim JCMT beam with a spa- 
tial shape spread over tens of beams, high signal-to-noise ratio SZ 
measurements of targeted clusters should be possible. Weaker sta- 
tistical measurements of less massive clusters may also be possible 
in 'blank sky' surveys. Despite the improvements in mapping struc- 
ture over many arc minutes, SCUBA-2's resolution will be the same 
as scuba's, and hence lensed sub-mm sources will still be an is- 
sue. However, with SCUBA-2's much higher sensitivity, their iden- 
tification and removal should be simpler. And obviously the avail- 
ability of deep data at other wavelengths will be crucial for extract- 
ing lensed sources. Candidate lensed sources in cluster fields can 
be followed up with pointed observations with much higher reso- 
lution using existing interferometers, and eventually the full power 
of ALMA. When point source free high frequency SCUBA-2 data 
are combined with those from the next generation of dedicated SZ 
effect experiments very powerful constraints on the full SZ effect 
spectrum in hundreds of clusters should become commonplace. 



4.3 Future experiments 

It is problematic for this type of measurement that, for a given clus- 
ter, background source lensing and the amplitude of the SZ effect 
are highly correlated. In other words, a massive cluster with a large 
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Figure 6. Comparison of the measurements in tliis work with those available in the literature. Each panel shows the results for one cluster; the y axis gives 
the measured yo and its associated 68 per cent confiden ce interval as Quoted in the refe rences as labeled o n the x axis. Thes e references are as fol lows: 
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the weighted mean of the decrement measurements (i.e. triangles) only. Interestingly, the points measured using the SZ increment tend to lie at larger i^o than 
their decrement counterparts; this may be indicative of point source contamination. 
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APPENDIX A 

Table 4: Archival cluster data sets. Listed are the clusters, JCMT project IDs, 
dates and integration times. We also give the mean and scatter of the optical 
depth for the uncorrupted parts of the data sets. For sets with insufficient Tgso 
measurements to measure the scatter, the ^(rgso) column is left blank. 



Cluster 


Project IDs 


Observation dates 


Integration 


(t850> 


'7(t"850) 








time (ks) 






Abel! (S8 


MOI AT ini 


Feb 77 7001 and Tan 4* 10 900? 

± tU. Z-V/v/ -L tlllU Jail. ^ OC A v/} 


12.2 


0.40 


0.04 


Abell 85 


M98AU60 


Apr. 19, Jun. 07, 1998 


9.7 


0.46 


0.08 




M98SU60 


Jun. 16, 1998 


1.0 


0.29 






Combined 




10.7 


0.44 


0.09 


Abell 209 


MOlAUOl 


Aug. 22, 23 & Nov 23, 2001 


24.3 


0.43 


0.05 


Abell 222 


M98AU17 


M. 02 & Nov. 29, 1998 


9.6 


0.20 


0.12 


Abell J /U 


U45 


T,,l AO f\A Q AC 1 AAT 

Jul. 02, 03, 04 & 05, 1997 


25.7 


0.26 


0.06 




JVl9/DJrll3 


/- -\ , nn 1 AAT 

Oct. 30, 1997 


6.3 


0.27 


0.08 




WLy / DXJJJ 


AiifT HQ 1QQ7 


7 7 


n 77 


n 07 






AiifT OA 9S Ifx Rt 11 1QQQ 


6Q 7 


n 77 


1 7 






and Nov 12 & 13, 1999 










M00BH04 


Nov 27 & 29, 2000 


5.1 


0.36 


0.06 




MOOBUHFLEX 


Dec. 08 & 09, 2000 


21.8 


0.37 


0.11 






and Jan. 07 & 09, 2001 










Combined 




135.8 


0.26 


0.12 


Abell 383 


M00BU18 


Feb. 28, 2001 


3.8 


0.26 


0.02 




MOlAUOl 


Feb. 22, Mar. 02 and lul. 26, 2001 


7.0 


0.40 


0.01 




M02BH46A 


Ian. 06, 2003 


5.1 


0.112 


0.01 




M03AH15A 


Feb. 18 &20, May 19 


18.0 


0.17 


0.05 






and lun. 03, 2003 










M03BH25A 


Aug. 19, 2003 


2.6 


0.30 


0.01 




Combined 




36.6 


0.20 


0.09 


Abell 478 


M97AN15 


Sept. 07 and Dec. 06, 1997 


9.0 


0.71 


0.39 




M98BU67 


Mar. 10 & 11, 1998 


12.8 


0.22 


0.03 




MQQRN1 ^ 




7 (S 


0.27 


01 




Combined 




24.4 


0.49 


0.37 


Abell 496 


M98BN02 


Nov. 14, 24, & 25, 1998 


20.1 


0.22 


0.07 




M99BN13 


Aug. 07 & 08, 1999 


15.4 


0.28 


0.03 




Combined 




35.5 


0.25 


0.06 


Abell 520 


M98AU17 


lul. 24, 25, 31 & Nov. 29, 1998 


30.6 


0.27 


0.07 




M98BN02 


Ian. 30, 1999 


11.5 


0.26 


0.03 




M00BN12 


Dec. 19, 20 & 21, 2000 


19.2 


0.31 


0.04 




Combined 




61.3 


0.29 


0.06 


Abell 586 


M98AU17 


Apr. 21, Sept. 12, Oct. 31, 


21.1 


0.19 


0.08 






&Nov. 01, 1998 
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Nov. 17, 2000 


1.3 


0.31 






M00BU18 


Nov 17, 2000 


6.4 


0.32 


0.01 
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7.7 


0.32 


0.01 


Abell 773 


MOlAUOl 


Feb. 22 and Mar. 03, 2001 


23.3 


0.39 


0.03 


Abell 780 


M98BN02 


Sept. 20, 1998 


2.6 


0.31 




Abell 851 


M98AU39 


Mar. 12 & 13, 1998 


30.1 


0.23 


0.15 




M00BH05 


Dec. 27, 28, 29, 30 & 31, 2000 


39.7 


0.25 


0.05 






and Ian. 01 & 09, 2001 










MOIAUHFLEX 


Mar. 30 & 31, 2001 


10.3 


0.25 


0.05 




Combined 




80.1 


0.25 


0.08 


Abell 963 


M00AU09 


Feb. 28 and Mar. 1 & 2, 2001 


9.6 


0.33 


0.03 



Continued on next page 
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Table 4 - continues from previous page 



Cluster 


Project IDs 


Observation dates 


Integration 


\ ' Sou / 


^ \ ' Sou / 








time (ks) 








iVlUUlJ U 1 


INOV. 1 / , zuuu 


s 1 
J. 1 


'X 1 






MOlAUOl 


May 7, 2001 


5.1 


0.48 


0.07 




M02BH46A 


Dec. 6, 22, 28, 29 & 30, 2002 


57.3 


0.19 


0.06 






and Jan. 06, 2003 










M03BH25A 


Oct. 03, 2003 


2.6 


0.30 






M04AH19A 


Mar. 9, 2004 


3.1 


0.33 


0.01 




Combined 




82.8 


0.24 


0.10 


Abell 1689 


M98AU17 


Apr. 21, 1998 


8.2 


0.18 


0.02 




M98AU27 


Feb. 26, 27 & Nov. 25, 1998 


23.0 


0.14 


0.04 




M99AU46 


Dec. 27, 1999 


7.7 


0.27 


0.06 




M99BN13 


Dec. 06, 07, 08 & 14, 1999 


35.2 


0.37 


0.15 




SCUBA 


Jun. 30, 1999 


6.4 


0.22 


0.06 




Combined 




80.5 


0.26 


0.14 






iviay ij tx, iu, z.uui 


1 A. 8 




u.uu 


Abell 1835 


M97BI16 


Jan. 31 & Feb. 01, 1998 


13.4 


0.17 


0.04 




M98AU39 


Apr. 03, 1998 


9.6 


0.20 


0.02 




Combined 




23.0 


0.18 


0.04 


AhpH IQM 


MQR ATT97 


Fph 1(\ 1QQ8 


J.O 


n 1 9 

U. i z 


n ni 


Abell 2163 


M98AC31 


Jun. 28, 29 & 30, 1998 


7.2 


0.17 


0.03 




M98AU27 


Feb. 26 & 27, 1998 


3.8 


0.12 


0.03 




Combined 




11.0 


0.15 


0.04 


Abell 2204 


M98AU60 


Feb. 10, 1998 


5.2 


0.38 


0.01 


Abell 2218 


M98AN14 


Mar. 20 & 22, 1998 


19.8 


0.39 


0.03 




MOOBNU 


Jan. 27&29, 2001 


19.2 


0.20 


0.02 




M00BU38 


Aug. 09, 11, 12& 13,2000 


55.9 


0.23 


0.10 




M01BN21 


Jan. 05, 06 & 07, 2002 


28.8 


0.34 


0.09 




1 rM n ri 
\^(J1 IIUIIICU 




1 9^ 7 

i Z J. / 


n 97 

u.z / 


n in 


Abell 2219 


M98AC41 


Apr 05, 1998 


16.1 


0.17 


0.02 




M00BN12 


Jan. 26, 2001 


14.7 


0.16 


0.04 




Combined 




30.8 


0.17 


0.04 


Abell 2261 


M98BC32 


Jan. 30 & Feb. 26, 1999 


4.9 


0.14 




/\Deil Z,jyyJ 




AiitT in 11 19 ^A 91 .^r 99 


L 1 .yj 


n 98 

U.Zo 


07 
U.U / 






iiriH Fipr 1 Q 1 QQ7 
aiiuucc.i7, vyy 1 










M98AU39 


Apr. 03, 1998 


6.6 


0.15 


0.01 




M99BH04 


Sept. 02, 03, 04, 08 


37.6 


0.28 


0.14 






and Nov. 12 & 13, 1999 










MOOBUHFLEX 


Dec. 08, 2000 


2.6 


0.44 


0.03 




MOIAUHFLEX 


Jun. 17, 18, 19, 2001 


24.2 


0.44 


0.08 




Combined 




98.7 


0.33 


0.13 


Abell 2597 


M98BN02 


Nov. 14, 15, 24, & 25, 1998 


23.0 


0.23 


0.09 


CI 0016+16 


M98AC41 


Apr. 05, Sept. 02, 1998 


13.0 


0.24 


0.01 




M98BC32 


Sept. 03,1998 


5.8 


0.42 


0.01 




M99BC41 


Nov. 11, 1999 


7.2 


0.31 


0.03 




M00BU18 


Oct. 13, 2000 


2.6 


0.37 






M02BC37 


Oct. 03, 2002 


2.1 


0.29 


0.01 




Combined 




30.6 


0.31 


0.10 


CI 0023+0423 


M00BN17 


Oct. 25, 2000 


16.6 


0.33 


0.04 




MOlBUOl 


Sept. 26, 30 & Nov. 15, 2001 


20.1 


0.28 


0.07 




NOlBUOl 


Dec. 08, 2001 


11.5 


0.24 


0.05 




Combined 




48.3 


0.29 


0.07 


CI 0024+1652 


M97BU33 


Aug. 14 & Dec. 19, 1997 


14.3 


0.20 


0.06 



and Jan. 26, 1998 

Continued on next page 
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Cluster 


Project IDs 


Observation dates 


Integration 


{t850> 


O-(t850) 








time (ks) 








M02BH46A 


Aug. 08, 2002 


0.4 


0.82 






M03AH15A 


Jun. 03, 2003 


2.6 


0.24 


0.02 




Combined 




17.3 


0.26 


0.18 


CI 0055-2754 


MOlAUOl 


Jan. 10, 2002 


2.6 


0.33 




CI 0152-1357 


M03BC31 


Sept. 28, 29 & Oct. 17, 2003 


13.4 


0.23 


0.03 


CI 0303+1706 


MOlAUOl 


Aug. 11, 12 & Nov 23, 2001 


13.2 


0.49 


0.07 


CI 0848+4453 


M00AN22 


Feb. 26, 27 & 28, 2000 


29.2 


0.27 


0.03 


CI 1455+2232 


M97BC45 


Jan. 04, 1998 


15.3 


0.15 


0.02 


CI 1604+4304 


M00AN22 


Mar. 03, 04, 21 & 22, 2000 


41.6 


0.17 


0.07 




MDflRNn 

iVlUUDi>l i / 


WCL. .ij, A\J\J\J Oc jail. AO, A\J\J L 


M 1 


U. 1 7 


u.u / 


CI 2129+0005 


M98AU27 


Feb. 27, 1998 


2.6 


0.11 


0.06 


CI 2244-0221 


U45 


Jul. 02, 03, 04 & 05, 1997 


23.0 


0.25 


0.06 




M98AU39 


Apr 03, 1998 


2.6 


0.12 


0.00 




Combined 




25.6 


0.24 


0.07 


MS 0440+0204 


M97BU33 


Aug. 09, 10, 12, 13 and 


31.4 


0.35 


0.13 






Sept. 22 & Dec. 19, 1997 










M98AU39 


Mar. 12 & 13, 1998 


7.0 


0.27 


0.03 




Combined 




38.4 


0.32 


0.13 


MS 0451-0305 


M98AC41 


Sept. 02, 1998 


1.5 


0.43 






M98BC32 


Sept. 03, 1998 


14.1 


0.41 


0.01 




M99BC41 


Nov 10& 11, 1999 


14.7 


0.31 


0.03 




M00BU18 


Nov 19, 2001 


5.1 


0.37 


0.01 




Combined 




35.5 


0.36 


0.06 


MS iOo4— Oozi 


M9oAL. 3 1 


Jun. Z6, l^Vo 


2.6 


O./Z 






M98AC37 


May 14, 19Vo 


5.1 


0.33 








JNov. 14, 2y 6L 3U, lV9o 


23.0 


0.24 


0.12 






Nov. lU & 11, iyyy 


14.7 


0.31 


0.03 




A/TQQRM1 


T\f^f- OA 07 OS 
UcL. UD, U/, UO, 


1 7 Q 


U.J J 


n 1 1 






anH M 1 QQQ 
aiiu iH", lyyy 










IVIUZ AiN i Z 


A ,^r- in 1/1 1 ^ 
Apr. lU, 14, ij. 


4D. / 


n ^7 


U.U/ 






anH IfS 700? 












Dec 93 94 96 77 9007 


75.6 


0.27 


0.11 






and Jan. 28, Apr. 23, 24, 25, 












26, 27, & May 19, 2003 










Combined 




167.7 


0.29 


0.10 


PKS0745-19 


M97AN15 


Sept. 12, 1997 


6.4 


0.36 


0.04 




M00BU18 


Sept. 20, 2000 


2.6 


0.13 






Combined 




9.0 


0.32 


0.10 


RXJ 1347-1145 


M98AI18 


May 31 & Jun. 01, 1998 


17.0 


0.56 


0.07 




M99AU77 


Jul. 03 & 04, 1999 


6.1 


0.42 


0.24 




M00BN12 


Dec. 19,20,21, 2000 


17.9 


0.31 


0.05 






and Jan. 25, 2001 










SCUBA 


Dec. 21, 2000 


0.4 


0.33 


0.04 




Combined 




53.0 


0.35 


0.16 


Zwicky 3146 


m98ac41 


Apr 08 & 09, 1998 


7.2 


0.16 


0.04 




m98bc32 


Jan. 28, 1999 


5.1 


0.18 


0.04 




Combined 




12.3 


0.16 


0.03 
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APPENDIX B 

Table 5: Positions, fluxes at 850 and 450 /^m, errors and statistics for 
the candidate sources. The values are determined by subtracting the fiducial 
JCMT PSF from the map centred on the nominal point source's position and 
then summing the squared residuals for the entire image. Sources which may be 
images of the SZ effect are marked; these are not subtracted before the SZ effect 
fit discussed in Section[3] Because the 450 fim SCUBA array is smaller than the 
850 fim array, some sources at the edge of the 850 fim map do not appear in the 
450 fim map; these are donoted with the label 'Off Edge' in the 450 ^m flux 
column. 



Cluster Source 


Q2000 


X 

O2000 


QCn iirr, 
OOU /ilTl 


oou /im 


X 1 






(hh;mm:ss) 


(-Hdd:mm:ss) 


Flux (niJy) 


Q /NT 
o/iN 


Degrees of Freedom 


Flux (mjy) 


Abell 209-1 


1:31:51.3 


-13:37:45 


10.7 


3.4 


171 / 177 


65+273 


Abell 370-1 


2:39:51.9 


-01:35:55 


20.5 


19.9 


403 / 177 


-2+19 


Abell 370-2 


2:39:57.6 


-01:34:50 


6.8 


7.8 


186/ 177 


1+12 


Abell 370-3 


2:39:56.3 


-01:34:25 


6.2 


6.5 


230 / 177 


6+6 


Abell 370-4 


2:39:54.0 


-01:33:34 


4.0 


4.1 


176/ 177 


1+3 


Abell 370-5" 


2:39:52.9 


-01:35:04 


3.0 


4.1 


142/ 177 


1+5 


Abell 383- r 


2:48:02.9 


-03:31:47 


5.5 


4.2 


156/ 177 


-23±22 


Abell 478-1 


4:13:26.9 


+ 10:27:41 


30.0 


16.2 


283 / 177 


11+33 


Abell 478-2 


4:13:28.2 


+ 10:28:03 


10.8 


4.2 


302 / 177 


0+20 


Abell 478-3 


4:13:23.4 


+ 10:26:54 


7.4 


3.5 


151 / 177 


3+24 


Abell 496-1 


4:33:37.6 


-13:15:45 


8.9 


7.8 


179 / 177 


3+12 


Abell 496-2 


4:33:42.9 


-13:16:08 


11.3 


4.1 


177 / 132 


8+47 


Abell 520-1 


4:54:09.6 


+02:54:48 


3.7 


3.3 


166/177 


0+13 


Abell 520-2 


4:54:06.8 


+02:56:33 


8.2 


3.6 


163 / 113 


101 + 199 


Abell 520-3" 


4:54:09.5 


+02:55:11 


7.5 


6.0 


200 / 177 


-1+21 


Abell 586-1 


7:32:17.6 


+31:37:34 


5.4 


3.5 


208 / 177 


1 + 16 


Abell 586-2" 


7:32:19.5 


+31:38:26 


7.3 


3.3 


183 / 177 


6+11 


Abell 773-1 


9:17:57.2 


+51:42:30 


7.7 


3.8 


172/ 177 


-53+115 


Abell 780-1 


9:18:05.3 


-12:05:41 


69.3 


17.1 


185 / 135 


75+58 


Abell 851-1 


9:42:55.5 


+46:58:43 


12.9 


16.0 


289 / 177 


7+15 


Abell 851-2 


9:43:05.8 


+46:58:33 


4.5 


3.3 


186/ 177 


39+27 


Abell 851-3 


9:43:05.1 


+47:00:10 


8.7 


4.5 


161 / 161 


-32+140 


Abell 851-4 


9:42:53.1 


+46:59:49 


5.5 


3.4 


193 / 177 


-2+13 


Abell 851-5 


9:42:48.4 


+46:59:22 


7.5 


3.8 


17/ 134 


35+35 


Abell 963-1 


10:17:08.0 


+39:02:59 


4.7 


4.5 


212/177 


-1+15 


Abell 963-2 


10:17:09.1 


+39:03:37 


8.5 


5.1 


196/ 177 


23+41 


Abell 963-3 


10:16:57.2 


+39:02:49 


4.5 


3.9 


193 / 177 


-14+57 


Abell 963-4" 


10:17:08.2 


+39:02:25 


4.6 


3.7 


218 / 177 


3+17 


Abell 1689-1 


13:11:34.5 


-01:20:20 


7.2 


6.4 


190/ 177 


9+28 


Abell 1689-2" 


13:11:29.2 


-01:20:45 


6.8 


6.6 


387/ 177 


3+18 


Abell 1835-1 


14:01:04.7 


+02:52:28 


13.4 


9.6 


225 / 177 


6+8 


Abell 1835-2 


14:00:57.5 


+02:52:51 


16.0 


10.2 


209 / 176 


-2+13 


Abell 1835-3 


14:01:02.1 


+02:52:43 


4.5 


4.1 


160/177 


1+6 


Abell 2204-1 


16:32:44.6 


+05:34:54 


16.8 


4.8 


175 / 177 


5+72 


Abell 2218-1 


16:35:55.5 


+66:12:58 


6.2 


8.7 


384/ 177 


2+25 


Abell 2218-2 


16:35:54.3 


+66:12:30 


10.8 


23.5 


187/177 


14+17 


Abell 2218-3 


16:35:51.2 


+66:12:07 


4.8 


8.7 


401 / 177 


8+29 


Abell 2219-1 


16:40:20.4 


+46:44:01 


9.6 


5.3 


189/177 


-1+14 


Abell 2219-2 


16:40:23.9 


+46:44:02 


8.2 


3.2 


198/171 


14+37 


Abell 2390-1 


21:53:33.2 


+ 17:42:51 


8.7 


6.9 


156/177 


6+5 


Abell 2597-1 


23:25:19.8 


-12:07:27 


14.5 


10.5 


202 / 177 


-2+14 


Abell 2597-2 


23:25:23.0 


-12:07:52 


6.6 


3.4 


156/177 


2+10 


Abell 2597-3 


23:25:23.6 


-12:06:51 


6.7 


3.8 


190/171 


-3+29 


CI 0016+16-1 


00:18:29.3 


+16:27:28 


27.0 


4.1 


151/90 


Off Edge 


CI 0016+16-2" 


00:18:33.2 


+ 16:26:10 


5.5 


3.7 


199 / 177 


-2+16 


CI 0023+0423-1 


00:23:54.5 


+04:23:03 


5.9 


5.7 


155 / 177 


13+24 


CI 0023+0423-2 


00:23:53.5 


+04:21:57 


4.5 


3.9 


174/177 


-10+29 



Continued on next page 
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V^l L4.13 L^i kj WLil 




<^2000 


850 iim 


850 ;;m 




4.50 //m 




f ViItttitti 


( +Hrl'TTlTTl '^JQ^ 


PliiY fmTv^ 

± iLijV I iiij y 1 


S/N 




PliiY fniTv^ 

± iL4.^ I iiij y f 


CI 0024+1652-1 


00:26:34.0 


+17:08:33 


21.0 


9.0 


196 / 162 


8+19 


CI 0024+1552-2 


r\r\ '^ n ^7 c 

00:26:37.5 


1 1 ^-7 r\r\ A c 

+17:09:45 


6.4 


3.0 


172 / 177 


2±38 


CI 0848+4453-1 


8:48:39.4 


1 A A O A'~l 

+44:53:47 


6.0 


4.1 


177 / 177 


—7+29 


CI 0848+4453-2 


8:48:39.9 


1 A A C A ^ ^ 

+44:54:22 


10.8 


6.6 


1 ^ O / 1 

173 / 177 


10+52 


CI 0848+4453-3 


O /I O 1 1 T 

8:48:31.7 


+44:53:51 


6.7 


5.0 


177 / 177 


— 14+31 


CI 1455+2232-1 


14:57:16.2 


+22:19:45 


7.8 


3.9 


207 / 177 


4±13 


CI 1455+2232-2 


14:57:16.0 


+22:20:45 


7.1 


3.1 


n o / 1 

208 / 177 


5±19 


CI 1604+4304-1 


16:04:26.0 


1 A o no n 

+43:03:57 


4.9 


4.4 


163 / 177 


12+20 


CI 1604+4304-2 


16:04:22.7 


+43:05:11 


5.4 


5.2 


212 / 177 


8+17 


CI 1604+4304-3 


16:04:23.2 


+43:05:29 


6.2 


4.2 


o o n / 1 

230 / 177 


—2+18 


CI 1604+4304-4 


16:04:32.3 


1 A ^ f\ A ^ O 

+43:04:23 


12.0 


3.6 


124 / 111 


Off Edge 


CI 1604+4304-5 


16:04:30.1 


1 (< o n /I .it 

+43:04:41 


5.3 


5.0 


O 'O / 1 

232 / 177 


—3+17 


CI 1604+4304-6 


16:04:28.8 


+43:05:11 


3.5 


3.4 


nn / 1 

209 / 177 


—6+18 


CI 1604+4304-7 


16:04:29.2 


1 JO r\ C A A 

+43:05:44 


3.5 


3.4 


168 / 177 


— 10+29 


CI 1604+4304-8 


16:04:23.7 


1 A '~i f\ A A f\ 

+43:04:40 


8.2 


9.0 


212 / 177 


18+20 


CI 2244—0221-1 


^ ^ .i ^ too 

22:47:13.8 


r\£\ f\/' r\f\ 

—02:06:09 


8.8 


3.8 


o n ^ / 1 

204 / 177 


6+18 


CI 2244—0221-2 


22:47:16.3 


—02:05:48 


4.6 


3.3 


197 / 167 


2+21 


CI Zlii— 0221-3 


22:47:12.3 


riri r\c on 

—02:05:39 


6.4 


3.7 


206 / 177 


7+21 


CI 2244—0221-4 


22:47:12.3 


n n n /I ^ 

—02:04:25 


6.5 


3.1 


148 / 131 


Off Edge 


CI 2244—0221-5 


22:47:10.1 


no nc n 

—02:05:57 


9.0 


6.0 


174 / 177 


6+28 


MS 0440+0204-1 


4:43:05.4 


1 nn 1 n o c 

+02:10:35 


5.0 


3.1 


161 / 177 


— 1+68 


MS 0440+0204-2 


4:43:06.8 


+02:10:26 


6.0 


4.1 


1 nn / 1 

199 / 177 


1+38 


MS 0440+0204-3 


4:43:14.6 


1 nn nn cn 

+02:09:59 


6.8 


3.1 


1 1 t A A 

152 / 144 


5+64 


MS 0440+0204-4 


A A ^ no c 

4:43:08.5 


1 nn 1 n ot 

+02: 10:37 


4.2 


3.4 


1 1 It T~l 

211/ 177 


—6+49 


MS 0451—0305-1 


4:54:09.2 


— 03:01:56 


8.3 


3.4 


163 / 177 


38+96 


MS 0451—0305-2 


4:54:07.1 


n o nn /i ^ 

—03:00:42 


7.9 


3.1 


164 / 172 


— 1+94 


MS 0451—0305-3 


4:54:10.6 


n o n 1 1^ A 

—03:01:24 


10.6 


5.4 


■^nn / 1 

200 / 177 


— 12+56 


MS 0451—0305-4 


A C A 1 O 

4:54:12.8 


n o nn ^ 

—03:00:52 


13.8 


7.9 


249 / 177 


2+60 


MS 0451—0305-5 


4:54:15.6 


no nn -in 

—03:00:40 


9.7 


5.0 


193 / 177 


—20+73 


MS 1054-0321-1 


10:57:03.7 


-03:37:12 


5.9 


4.5 


190 / 177 


3+13 


MS 1054-0321-2 


10:56:56.0 


-03:36:34 


5.3 


4.1 


152/177 


15+26 


MS 1054-0321-3 


10:56:56.5 


-03:36:11 


6.2 


5.7 


163 / 177 


25+37 


MS 1054-0321-4 


10:57:02.0 


-03:36:06 


6.0 


4.4 


180/177 


15+21 


MS 1054-0321-5 


10:57:00.2 


-03:35:39 


5.0 


4.1 


201 / 177 


31+20 


RXJ 1347-1145-1 


13:47:27.6 


-11:45:54 


15.1 


5.1 


140 / 177 


125+34 


RXJ 1347-1145-2" 


13:47:31.3 


-11:44:57 


11.4 


4.7 


191 / 177 


10+32 



"These sources are possibly images of the SZ effect and we do not subtract them before doing the SZ effect 
fit. Because in an number of cases the apparent source is resolved, the in some of these sources is poor. 
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APPENDIX C 

This appendix contains notes on the characteristics of individual 
clusters, with emphasis on cm/mm/sub-mm and X-ray measure- 
ments of point sources and extended emission in each. 



Abell 68: ICoorav et al.l l ll998h has studied this moderate redshift 
cluster as part of a 28.5 GHz study of point sources in galaxy clus- 
ters. 

Abell 85: This is a low redshift cluster with no apparent sub- 
mm point sources or extended emission. Measurements of the SZ 
effect in this cluster around 30 GHz have bee n published previously 
dUdomprasert et al.ll2004lBaker et alj|l994l) . 

Abell 209: This cluster is relatively poorly studied, particularly 
at (sub-)mm wavelengths. 

Abell 222: This is a moderate redshift cluster with no apparent 
sub-mm emission, although the SCUBA map is quite shallow. Be- 
cause the isothermal /3 model for this cluster had to be derived from 
archival X-ray data, it has not been considered a viable SZ effect 
field and has been removed from the analysis discussed in Section 
I3.3l and later. 

Abell 370: This field cont ains the brightest so urce in this cat- 
alogue, discussed in detail in llvison et"al1 ( Il998|). Our results are 
consistent with th ose found by other groups Jsmail et al 
iBargeret all 19991) . As shown in lGrego et allJaOOOh. the SZ decre- 
ment s hows quite complex structure in this cluster. iLaRogue et al.l 
( l2006h provid es a newer measure ment of the SZ effect using 
OVRO/BIMA. ICoorav et aP ( Il998h provides a point source list at 
28.5 GHz for this clu ster, and ISOCAM re sults at 14 and 7 jim have 
also been published jMetcalfe et al.ll2003h . 



Abell 383: ICoorav etal.l ( ll998h provides a 28.5 GHz point source 
list for this moderate redshift cluster. 

Abell 478: This cluster has b e en a popul ar target for SZ ef - 
fect experiments; iRadford et alj ( Il986h and IChase et"al ] l ll987h 
discuss early attempts to measure t he SZ effec t at 3 and 
1.2 mm, respectively. More rec ently, iMvers et al.l ( 1 19971) and 
iMason. Mvers & Readhea j (12001 ) discuss measurement s of the SZ 
effect using OVRO at 32 GHz. lUdomprasert et alj(l2004h 



Lancaster et al.l feOOSi) 



measures 
uses 



this cluster's SZ effect with CBI, and 
the VSA at 34 GHz as well. 

Unfortunately, this cluster is a poor target for SZ effects 
at sc uba's shorter wavelengths. iKnudsen. van der Werf & Jaff3 
(l2003h have presented the data discussed here and associate the 
bright point source near the centre of this field with a high redshift, 
dusty quasar. Previously, £ox, Bregman & Schombert ( 1995) dis- 
cuss IRAS measurements of this clust er at 60 and 1 /jm. 

IPointecouteau et al.l ( l2004h and Ide Plaa et al] ( |2004 discuss 
the X-ray properties of this cluster using recent XMM-Newton data, 
an d ISunetalJjiooal) presents similar measurements with Chandra. 

A bell 496: This is a low redshift cluster, discussed extensively 
in lDurret et al ]|2000l We dete ct a bright source as sociated with the 
central galaxy in the cluster ^Griffith et alJll994h . This cluster is 
therefore contaminated from the point of view of SZ effect mea- 
surement with the JCMT. 

Abell 520: IChapman et aP ( l2002h provide a previous analysis of 



the sub-m m data for this cluste r. The SZ effect has been measured a t 
30 GHz bvlReese et aljj200 A . with SuZIE II llBenson et al.l2004l) . 
and lCoorav et alj 1 19981) provide a list of point sources at 28.5 GHz 
in this field. 

Abell 586: A deep X MM-Newton observation is presented in 
iDe Filippis et al ] ll2003h . and 28.5 GHz d ata focusing on poin t 
sources are given in 



Coorav et all ( Il998l) . lUaRoque et all ( l2006h 



presents an SZ effect measurement at 30 GHz using OVRO/BIMA. 

Abell 665: This moderate re d shift c luster is another p opu- 
lar SZ effect target. iGrego et all ( 1200 ih , iReese et all ilOO^ and 
lLaRoQueet"aD(l2006h measure the SZ effect decrement at 30 GHz, 
and ICooravetalJ lll998l) measur e point sources at the same fre- 
quency. Also, Desert et al] l ll998l) discusses DIABLO SZ observa- 
tions at 2.1 and 1.2 mm. 

Abell 773: The clu ster has a number of SZ effect measurements. 
iBenson et alj j2004h present measurement of the SZ e ffect using 
SuZIE II. Both 'Grai nge et all l ll993l) . ISaunders et all ( l2003b and 
llones et al. (.2005) discuss a pro gram of SZ effect measu r ement s 



using the Ryle te lesco pe. Also. [ Carlstrom, Jov & Gregd ( Il996h . 
iReese et al. I J2002h and lLaRogue et al.l ilOO^ SZ effect measure- 
ments at 30 GHz in this cluster. Recentlv. lGovoni et aD (l2004l) pre- 
sented X-ray data from Chandra for this cluster. 

Abell 780: Abell 780 contains Hydra A, an extremely brigh t ra- 
dio source (AGN), in its central regions dBennett et al.ll2003l) . As 
Hydra A has a flux of about O.lJy at 850 fj,m, it is not possible to 
isolate the SZ effect in this cluster. 

Abell 851: Previous sub-mm meas urements of this clu ster are 
presented in lSmail et al.l | |2002|) . while lCoorav et al](ll998h present 
measurements of point sources at 28.5 GHz. iDe Filitjpis et all 
ll2003h find a complex X-ray structure in this cluster which exhibits 
two cores, each having a temperature of about 4.7 keV. We there- 
fore expect that the isothermal /3 model is a poor description of the 
shape of the SZ effect emission in this cluster, and exclude it from 
the analysis discussed in Section [33] 



Abell 963: This cluster is very well studied at optical wave- 
lengths, and as a part of a number of X-ray surveys. 

Abell 1689: Our map of Abell 1689 exhibits a great deal of 
extended sub-mm emission around the expected centriod of SZ 
effect emission, which is best explained as a sub-mm image of 
the cluster's SZ effect. This cluster has been a ver y popular tar- 



get fo r SZ effect measurements, including SuZIE II (iBenson et al 



20041) and OVRO/BIMA at the SZ decrement llHolzapfel et al 



1997*. 'Grego et al.' 200 iL iReese et al.|[200l LaRoque et aljbooji 
Coorav et al. ( 1998) discuss point sources at 30 GHz. 



Abell 1763: ICoorav et al.l d 19981) present point source measure- 
ments at 28.5 GHz for this cluster. Abell 1763 has been detected in 
a number of X-ray surveys. 

Abell 1 835: This field contains two known lensed sub-mm 
sources ( llvison etal]|2000l) and has also been imaged at 1.1 mm 
with BOLOCAM (D. Haig, private communication). In addition, a 
small amount of emi ssion in the centr al regions is detected at rea- 
sonable significance. lEdge etalJ ( ll999l) associate this flux with dust 
emission from the central galaxy; this is supported by the detection 
of CO emission in this galaxy. 
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and lCoorav et alj 



Grego et alj ( 2001 ) detect an SZ effect decrement at 30 GHz, 
1998 ) discuss p oint s ources detecte d at the same 



frequency. iMauskopf et al] ilOOd) and ISenson et all l l2003h both 
find data consistent with a sub-mm SZ increment in this clus- 
ter using SuZlE I/ll. Bo th Schmidt. All en & Fabiaiil l l200lh and 
iMai erowicz. Neumann & R eiprich ( 2002) present X-ray measure- 
ments of this cluster, finding a central ICM temperature of 4.0 and 
4.4 keV, respectively. Abell 1835 has a strong cooling flow, and 
thus exhibits a temperature gradient between the centre and the hot- 
ter outer regions. 

Abell 1914: This field only has about an ho ur of SCUBA inte - 
grat ion time, resulting in a very shallow map. iGrego et"al] ( 1200 ih 
and iLaRooue et alj 1 200d) give prior SZ effect measurements at 

discuss point sources in this clus- 



30 GHz, and 



Coorav et al 



ter. Also. ljones et alj l l2005h use the Ryle telescope to measure the 



SZ effect in this cluster. 



Abell 2163: This cluster is very well studied at mm wave- 
lengths, particularly in terms of the SZ effect. SuZIE I/II has 
measured this cluster a number of time s ([Wilbanks et al. [ 1 1994 
Holzapfel etaLlll997L iBenson et al]|2004h . iReese et alj ( l2002h and 



LaRogue et alj 1 200( ) present SZ effect measurements at 30 GHz 



using BIMA/OVRO. The DIABOLO experiment has been used to 
measure the SZ ef fect in this cluster at millimetric wavelengths 
jPesertetalJl 19981) . and PRONAO S measured the sub-m m SZ ef- 
fect for the first time in Abell 2163 jLamarre et"al]| 19981) . 

Sadly, the SCUBA data for this cluster are pathologically 
noisy; the map shown in Fig|2]has structure that is not related to as- 
tronomical emission. We have attempted to understand this noise, 
but no obvious cause has been found. We therefore present the 
850 /im map, but do not analyze the data beyond that stage. 

ICoorav et a l.' fl99^ p resent 28.5 GHz point sources in this 
cluster, and Icovoni et afl (l2004h presents recent data from the 
Chandra observatory. Chapman et alj ( |2002|) have discussed some 
of the SCUBA data presented in this work. 

Abell 2204: Th e SZ effec t in t hi s cluster has previ ously 
been measured bv [ Holzapfell (Il996l). iBenson et alj ( |2004 and 
lUaRogue et alj ( |2006|) . while Coorav et alj 1 19981) give a discussion 
of 28.5 GHz point sources in this cluster. 

Abell 2218: This is among the best studied clusters in the sky, 
and is particularly popular for studies of gravitational tensing. In 
fact, one of the candidates f or the highest redsh ift source (z ~ 7) 
known is in this cluster field jKneib et al.ll2004al) . 

iKneib et alj l l2004bh present the discovery of SMM 
J16359+6612, a high redshift, lensed sub-mm galaxy. This 
source is present in our maps, but because iRneib et alj ( l2004bl) use 
a smaller pixel size, more structure is evident in their maps. This 
cluster is not used in our final SZ sample because the lensed source 
is too complex to model and remove from the data. 

At wavelengths close to SCUBA's. fsheth et alj j2004h discuss 
the de tection of CO from SMM J16359+6612, and iBiviano et alj 
( l2004h present an ISOCA M study of star forming galaxies in Abell 
2218. Egami et"al] 1 20051) prese nts Spitzer observations of this clus- 
ter as well. lCoorav et al. list radio point sources at 28.5 GHz 
in this cluster. 

Abell 2218 was also a popula r target for early sea rches for 
the SZ effect like that present ed in iRadford et aO ( Il 986i) and cul- 
minating in work like that in Ijones et al. ( 1 1993 ). More recently, 
this cluster has been imaged by the Effelsberg 100 m telescope at 



10 GHz dUvaniker et aO 19971) , the Nobeyama telescop e at 36 GHz 
jTsuboi et alj Il998l). SuZIE II l lBenson et all l2004t) . the Ryle 
telescope jjones et alj |2005[). and OVRO/BIMA around 30 GHz 
jGrego et al.l200lLlReese et alj200irLaRooue et al.l2006l) . It is not 
obvious that this cluster is a particularly good target for SZ obser- 
vations above these frequencies due to SMM J163 59+6612. 

Recent X-ray studies of this cluster include Machacek et al] 

l l2002h using Chandra and IPratt. Bohringer & Finoguenovl ( |2005|) 
using XMM-Newton. 



Abell 2219: IChapman et alj (l2002h provide a previous analysis 
of these sub-mm data, while lCoorav et al. Il ll998h discuss 28.5 GHz 
point sources in this cluster. 

Abell 2261: The SCUBA integration time on this cluster is 
only a bout 1.5 hours, resulting in a shallow map. IChapman et al] 
l l2002h have presented a previous analysis of this sub- mm data. The 
SZ effect has been mea sured both at the dec rement ( Grego et alj 
2000llReese et aljl2002l) and around the null ( IBenson et alj|2003l 
Benson et alJl2004l) . ICoorav et alj ( ll998l) discuss point source con- 
tamination near 30 GHz. 



Abell 2390: These sub-mm data are discussed i n iBarger et al] 
( Il999h . and lEdge et all ( Il999h . while IBenson et al J ( |2003|) provide 
a measurement of the SZ effect in this cluster using SuZIE II. 

Abell 2597: Abe ll 2597 contains a bright AGN 

dVeron-Cettv & Verod 1200 ih which is imaged in this sub- 
mm data set. It is therefore no t suita ble for SCUBA SZ effect 
observations. lUdomprasert et all ( l2004h present a SZ effect study 
of this cluster using CBI. 



CI 0016+ 16: This cluster is probably the most popu- 
lar SZ eff ect target in the sky; in addition to our previ- 
ous work dZemcov et al] l2003ln . a number of experiments 
have measured the SZ eff ect in this cluster. Early mea- 
sure ments are discussed bv iBirkinshaw^ Gull & Moffet 1 1981 ) 



and iBirkinshaw. Gull & Hardebeckl l ll984i) . IBenson et alj I2OO3I) 
and Benson et alj 1 20041) discuss measurements using SuZIE 
II on the Caltech Sub-mm Observatory. Measurements us- 
ing the OVRO/BIMA system around 30 GHz are discussed 



ICarlstrom. Jov & Gregc (Il996l) , [Hughes & Birkinshawl 
iReeseetal.' 1200(f). iGrego et alj ll200ll). iReese etal. " 



1998), 



2002), 



.LaRoque et al. (2003il. and iLaRoque et all l l2006h . IPesert et alj 



1 1199 8) pre sent SZ effect mea surements using DIABOLO on 
IRAM, and lGrainge et alj ( l2002h used the Ryle telescope for their 

SZ measurements. 

IChapman et alj (l2002h alrea dy presented an analysis of these 
sub-mm data, while ICoorav et al . (1998) discusses measurements 
of radio point sources at 28.5 GHz in this cluster. This clus- 
ter has also been targeted w ith a number of X-ray telescopes; 
IWorrall & Birkinshawl f2003h presents recent measurements usmg 
XMM-Newton. 

CI 0023 + 04 23: This high z cluster was origirially di scovered in 
the survey of ICunn. Hoessel & Okel ( ll986l) . lBesll ( I2OO2I) has previ- 
ously analyzed these sub-mm data. As no isothermal /3 model pa- 
rameters are available or derivable from archival data, an SZ effect 
measurement cannot be performed with these data. 



Th ere are approximately 3 times more data in the work presented here 
than in IZemcov et all <2003ll . 
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Cl002i + 1652 
cussed 



i + lbt)Z: These S CU BA data have previo usly been dis- 
in lSmail et all ( l2002h and lFraver etal] hoOOh 



CI 0055 — 2754 ■ This is a relatively poorly studied cluster; due 
to both a lack of a well-defined SZ model, and the short integration 
time on this field, it has been rejected from our SZ sample. The map 
is presented in Fig.|2] 



CI 0152-1357: iMaughan et"al] ( |2003|) discuss X-ray imaging 
of this cluster, and show that it is composed of two clumps. Our 
SCUBA data image d the northern cl ump, which has an X-ray tem- 
perature of 5.5 keV. IJoy et al.l ( 1200 ll) have previously measured the 
SZ effect in this cluster at 30 GHz. 

CI 0303+1706: This cluster is relatively poorly studied, partic- 
ularly at mm wavelengths. 

CI 0848+44 53: An an alysis of these data has been presented 
previously bv lBestlJiool . 



CI 1455+2232: IChapman et al.l (l2002h presen t a pre vious anal- 
ysis of these sub-mm data, while ICoorav et al.] ( Il998h discuss a 
28.5 GHz survey for point sources. 



CI 1 604 +4304: iGunn. Hoessel & Okd l ll986l) firs t detected a 
set of galaxy clusters in this region of the sky, and iLubin et^ 
( l2000h showed that this c l uster forms part of a larger supercluster 
with CI 1604+4321. lBestl ( l2002l) previously presented this sub-mm 
data. As the isothermal 13 model had to be estimated in this cluster, 
it has not been considered in the analysis presented in Section [33] 



CI 1604 +4321 : This is the partner of CI 1604+4304 and the 
same comments apply. 

a 2129+0005 : ^This cluster hosts a bright AGN 

( ICrawfordet"aLlll999h . and so is not suitable for SCUBA SZ 
effect measurements. 



CI 2244 - 0221 : ISmail et alj ( 1 19971) present an analysis of the 
sub-mm emissi on in this cluster based on these SCUBA data. 
lota et al] ( Il998h surveyed this region with ASCA, and find a best 
fitting isothermal model with /3 — 0.30. This value of /3 is phys- 
ically acceptable for the X-ray emission, but produces a divergent 
model for the SZ effect. We therefore use the classic King profile 
(13 = 2/3) for this cluster. However, this model has a very large 
core radius, and poorly describes these sub-mm data. As the model 
is such a poor description of the known properties of this cluster, 
we have not considered it in the analysis discussed in Section [33] 
and later. 

MS 044 + 0204 ■ These S CU BA data have previo usly been pre- 
sented in lBarger et al] ( Il999h and lSmail etall J2002h . 

MS 0451- 305: This is a high redshift cluster in which 
iBorvs et al.l ( |2004|) have found a resolved, bright submm source. 
The sa me group presented a map made with some of these data 
earlier dChapman et"ai]|2002h . but did not have low enough noise 
to identify the source. While this source is interesting in its own 
right, it precludes the type of SZ analysis discussed in this paper. 

This cluster is another popular target for SZ effect observa - 
tions. It was discussed in our earlier work IZemcov et ah I i l2003h . 
but was rej ected from th a t SZ sampl e due to the l ensed sub- 
mm source. iBenson et al.l l l2003h and iBenson et al.l l l2004h dis- 



cuss observations of this cluster with SuZIE II. MS 0451—0305 
is also a popular target f or OVRO/BIMA at 30 GHz jReese et al.l 
20001 iGrego etall 1200 ll iReese etal.ll2002l lUaRogue et all l2003l 



LaRooue et al.l2006l). 

ICoorav etal .' ('1998*) d escribe a sea rch for radio point sources 
at 28.5 GHz, and Donahue et al. I l l2003l) presents recent Chandra 
maps for this cluster. 

MS 1054 — 03: This is a high redshift, hot cluster of galax- 
ies. It has been targeted for SZ effect observations by 
SCUBA jZemcov et al]|2003l) SuZIE II llBenson et al.ll2004h and 
OVRO/BIMA jjoY et al.l200lllLaRogue et al.l2006l) . It is also well 
studied in the X-ray, mos t recently with C handra l ljeltema et al.l 

l200lh and XMM-Newton jGioia et alj2004h. 

Both ICoorav et ah! ( 1 19981) and Best et al.' ('2002') provide ra- 
dio source lists for MS 1054— 03.iChapm an et al.. (2002i) previously 
discussed some of the data in the set presented here. A large frac- 
tion of t he total integration tim e came later, and those data are pre- 
sented in lKnudsen et al.l ( l2005h : they present a deep source list, but 

not the map itself. 

Some of the data for this cluster also appeared in lZemcov et al.l 

( l2003h : the best fit SZ effect given there is consistent with that given 
in this work. 



PKS 0745-19: lEdgd ( 1200 ih has presented this data as part of a 
progra mme to measure CO emission in galaxy clusters. [Hicks et al] 
( I2OO2I) presents recent Chandra data for this cluster. 



RXJ 1347-114 5: RXJ1347- 



1145 is t he brig htest cluster in the 
iSilOOl!)] 



ROSAT sample: ! Allen. Schmidt & FabianI (12002!) present measure- 
ment of this cluster with Chandra. They find a central X-ray temper- 
ature of 7.0 keV that rises to 16 keV away from the cluster centre. 

This cluster has been a target for ma ny SZ effect 
exper i ments from 30 G Hz to 350 GHz dKomatsu et all 



im 



Pointecouteau et al.l 



199^ 



Pointecouteau et al.l 



Komatsu et al.l I2OOIL iReese et al] |2002| . IBenson et al.l 



2001 



2004 



LaRoque et al I I2OO6I) . It exhibits what appears to be among the 



only images of the SZ effect at 350 GHz in our sample , although 
we find a different flux than that given bv lKomatsu et aLl ( [l999D . 

RXJ 1347—1145 exhibits the only detected 450 ^m source in 
this catalogue. Unfortunately, it is not clear what source this corre- 
sponds to at other wavelengths. Due to its proximity to the centre 
of the cluster, this source may be confused in SZ effect experiments 
with larger beams than that of SCUBA on the JCMT. 

Zwicky 31 46: These data have previously been analyzed and dis- 
cussed bv lChapman et alj(l2002l) . The SZ effect has been measured 
by SuZIE I/II (lHolzapfel"l996'. 'Bens on e7ai]|2003l . IBenson et al.l 
[2004), while Cooray et al. C199&) discuss radio point sources at 
28.5 GHz in this cluster. 

This paper has been produced using the Royal Astronomical Soci- 
ety/Blackwell Science KTgX style file. 
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